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PREDICTING SEEPAGE UNDER DAMS ON MULTI-LAYERED FOUNDATIONS 


Paul H. Shea! and Harry E. Whitsett? 
(Proc. Paper 1727) 


SYNOPSIS 


During the design of the Central and Southern Florida Project it became 
necessary to develop methods for predicting rates of seepage under dams and 
levees on multi-layered foundations since there were no published mathemati- 
cal or graphical methods which were applicable. The methods developed are 
presented in this paper. A description of the pumping test procedure used to 
implement those methods is included. 


INTRODUCTION 


Determination of the quantity and effect of seepage under dams and levees 
on pervious foundations may be accomplished quite simply when the foundation 
contains a single pervious layer which is open to the reservoir or river. 
Where the pervious foundation is blanketed by a relatively impervious layer, 
underseepage may be computed by the method described by Bennett? or may 
be derived from a flow net when differences in permeability are not too great. 
In cases where more than one pervious layer and more than one impervious 
layer is present, the construction of a flow net is so time-consuming that it 
becomes impractical. Even in cases where there is only one pervious and one 
impervious layer, a flow net is impractical if there is a very great difference 
in permeability between the layers because of the need for drawing either 
Note: Discussion open until January 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1727 is 


part of the copyrighted Journal of the Soil Mechanics and Foundations Division, 


Proceedings of the American Society of Civil Engineers, Vol. 84, No. SM 3, August, 
1958. 
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extremely large or minutely small “squares”. Solution of such problems by 
mathematical processes was found to be possible, however, provided that the 
field conditions conform reasonably well to the following assumptions: 


(1) Each blanket and each aquifer is a homogeneous medium of uniform 
thickness. 

(2) Each blanket has a constant vertical permeability and each aquifer has 
a constant horizontal permeability. 

(3) Flow is laminar and therefore conforms to Darcy’s law. 

(4) Flow is vertical in the blankets. 

(5) Head loss accompanying vertical flow in the aquifers is negligible as 
compared to vertical head loss in the blankets. 


If the results of borings indicate that the geologic conditions at a given site 
conform generally to assumption (1) above, conformity to the other as- 
sumptions may be checked by one or more field pumping tests. Such tests are 
the only practicable means for determining the transmissibilities of the vari- 
ous strata. Laboratory testing of so-called “undisturbed” samples would pro- 
vide a very unreliable basis for the required analysis. Aside from the 
enormous number of samples which would be required to obtain representa- 
tive permeabilities, there is the added difficulty that those materials which 
are most pervious are the most difficult to sample and also the most difficult 
to prepare for testing. In the case of rock foundations, where leakage occurs 
along open joints and bedding planes and through solution holes, laboratory 
testing of rock samples would be futile. 

If foundation conditions are so non-uniform that the observed radial flow 
pattern to a pumped well cannot be analyzed, there is no method by which flow 
through such a foundation under prototype conditions can be determined. On 
the other hand, if the radial flow pattern can be analyzed, the flow under a 
dam or levee on the same or similar foundation can be predicted with reason- 
able accuracy. 

From the foregoing discussion it can be seen that there are three pre- 
requisites to the computation of seepage quantities under water-retaining 
structures on multi-layered foundations. They are (1) direct observation of 
radial flow in the foundation by means of a pumping test; (2) a method for 
analysis of radial flow in a multi-layered foundation; and (3) a method for 
analysis of flow under a dam or levee on the same foundation. 


Development and Use of Procedure for Underseepage Estimates 


The necessity for development of a method for predicting seepage quanti- 
ties under levees on multi-layered foundations arose during the design of the 
Central and Southern Florida Project. That project (still incomplete) will 
eventually be an integrated system of water-control works for flood control, 
drainage, and water-conservation—covering all eastern Florida from the lati- 
tude of Orlando southward. Included in the project will be hundreds of miles 
of canals and levees, numerous small dams and spillways, and several very 
large pumping stations. 

The project provides for the retention of excess rainfall in conservation 
areas to be ringed with levees. Water thus stored could be used, when needed, 
for irrigation. Obviously, it would be desirable to keep seepage losses at a 
minimum from those areas. Unfortunately, much of the reservoir area is 
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underlaid by very porous and permeable limestones extending from—or near— 
the ground surface to depths of 50 feet or more. The limestones are perfo- 
rated with solution channels ranging from a fraction of an inch to 6 in. or 8 in. 
in diameter. 

Throughout most of the area, the uppermost limestone is much less pervi- 
ous than the remainder because of the sealing of pores and solution channels 
by peat, marl, or redeposited calcite. The thickness of the relatively tight 
layer varies from a few inches to 2 feet or 3 feet. A further complication is 
the presence of horizontal beds of relatively impervious fresh-water lime- 
stone dividing the more pervious marine limestone into two or three aquifers. 

Prevention or reduction of underseepage by grouting the porous rock or 
providing cutoff walls is not economically feasible. The cost of such 
treatment would be many times the cost of the levees. Canals, pumping 
stations, and control structures in the areas that will be affected by the 
seepage must be designed with capacities that will include underseepage flow. 
A reasonably accurate estimate of the quantity of underseepage that will occur 
when water is stored in the conservation areas is required for the design of 
these structures. 

Soon after work on the project began it became apparent that none of the 
published methods for the determination of underseepage was applicable to 
the multi-layered condition encountered. The development of the computation 
methods described in this paper was necessary to the design of the project. 

Guidance in the early stages of the underseepage investigations was 
furnished by P. T. Bennett4 and R. A. Barron.5 The pumping test procedure 
originally set up was followed with only minor modifications throughout the 
investigation. The derivations and solutions of the differential equations for 
flow in systems having more than one aquiclude and aquifer were done by the 
junior author and are included as appendices. Those for flow in a single- 
blanketed aquifer were published or implied in Bennett’s paper referenced 
above. 

The estimates of underseepage used in design were based on pumping tests 
performed at rather widely spaced locations along the levee alinements. 
Interpolation between tests locations was based on geologic studies of core 
samples and the results of “recharge” tests made by pumping into core bor- 
ings along the alinements. 

Predicted underseepage quantities have been checked by metering the 
continuous borrow pits on the landside of several of the levees. The metering 
results have been erratic but they have indicated that seepage quantities esti- 
mated from pumping tests have good probable accuracy. Metered quantities 
have varied from 60 to 110 per cent of those estimated. Comparisons of com- 
puted and measured underseepage are shown on Figs. 1 and 2, “Measured 
underseepage” is the increase in flow between two points in the continuous 
landside borrow pit which parallels each levee. That flow does not necessari- 
ly represent all the underseepage since some of the seepage may by-pass the 
borrow pit and particularly at low stages, some flow may be contributed by 
depletion of ground water storage. Fig. 1 shows very good agreement between 
the computed and measured values. In Fig. 2 the metering results are very 
erratic. The apparent average of the measurements is about 70 per cent of 
the computed underseepage. After several years of metering, the seepage 


4. Engineer, Missouri River Division, C. of E. 
5. Engineer, Office, Chief of Engineers. 
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quantities computed from the pumping tests are still considered the most ac- 
curate data for use in hydrologic studies. 


Details of Procedure 
General 


As stated above, the process of making an underseepage estimate consists 
of three steps: Performing a pumping test or tests; computing transmissibili- 
ties from the results of the test; and computing prototype underseepage from 


the transmissibilities. Uplift pressures at the base of the top blanket or aqui- 
clude can also be computed. 


The Pumping Test 


The proper performance of the pumping test is of paramount importance. 
Adequate observational data are essential if a reasonably accurate prediction 
of underseepage is to be made. A procedure which has been used successfully 
by the U. S. Army Engineer District, Jacksonville, Florida, is described in 
some detail below. 

Preliminary Exploration.—Small diameter borings are put down to de- 
termine the ceneral geology of the area. The depth to which pervious materi- 
als extend and the elevations of impervious layers are determined. Sites for 
one or more pumping tests are selected at points where conditions are best 
suited for the tests. If possible, tests are performed at locations where the 
effects of streams, drainage from hillsides, and other features which might 
affect the results are unimportant. 

The Pump Well.—The pump well is drilled to the full depth of the pervious 
formations. It is screened and gravel-packed through unconsolidated sedi- 
ments and uncased through rock except that slotted casing is installed in zones 
where caving is anticipated. A clay or grout seal is placed at the surface to 
prevent entrance of surface water. 

Placement of Piezometers.—Depths of observation wells or piezometers 
are based on a study of the samples recovered in the preliminary borings and 
on a vertical velocity-meter traverse made up and down the well during pump- 
ing to determine the distribution of inflow to the well. A set of piezometers 
is installed in each layer which contributes an appreciable portion of the in- 
flow. Piezometers are located along straight lines passing through the pump 
well. Spacing of piezometers is selected to give equal increments of the log 
of the radial distance from the pump well. The length of the lines must be 
great enough to define both the straight and curved portions of the curve of 
drawdown versus the log of the radial distance. An estimate of that length is 
made for the initial installation and the lines are extended, if necessary, after 
the drawdown curves are plotted for the first pumping of the well. Standard 
well points are used for piezometers placed in sand or gravel beds. For 
piezometers in rock, open end pipes are used. The bottom of the pipe is 
placed one foot above the bottom of a drilled hole and is grouted into place. 
All piezometer pipes extend far enough above the ground to prevent debris 
and water from entering. Piezometers are flushed with clear water after 
installation to make certain that they are clean and open. 

Methods of Measuring Drawdown.—In very pervious formations, drawdowns 
are small and very accurate methods of water-level determination are re- 
quired. If the water level in the piezometers is high enough, a point gage 
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connected in series with a small battery and voltmeter gives excellent results. 
When the point of the gage makes contact with the water surface, completion 
of the circuit is indicated by the voltmeter. Readings accurate to .001 foot 
can be obtained by using the verneir on the point gage. When the water level 
is too low to permit use of the point gage, another electrical water-level indi- 
cator can be used. This instrument makes use of two insulated wires, with 
tape markings for measurement of depth, connected in parallel with a battery 
and neon lamp. When the gage comes in contact with the water, the circuit is 
shorted out, extinguishing the lamp. Readings accurate to .01 foot can be 
made with this instrument. 

Test Procedure.—Before pumping the well, the distance from the top of the 
pipe to the water surface is measured in each piezometer. The static levels 
thus measured furnish the datum from which drawdowns observed during 
pumping are plotted. True elevations are not required. After static levels 
are established, the time lag for each piezometer is checked by filling the 
pipe with water and observing the time interval required for the water to re- 
turn to the static level. The time lag is nearly always negligible. After static 
levels and time lags are known, the pump is started, and pumping is continued 
at a constant discharge rate. Piezometers are then read at regular intervals 
until little or no change in drawdown can be observed between successive 
readings. Drawdown curves are plotted as soon as the first pumping test is 
completed. An examination of the curves indicates where additional piezome- 
ters, if any, are needed, and whether the pumping rate is in the optimum 
range. If the pumping rate is too high, drawdowns in the piezometers nearest 
the well are excessive, indicating turbulent flow. If it is too low, drawdowns 
in the piezometers farthest from the well are too small to plot accurately. 
Any necessary chaages in piezometer lay-out are made before the well is 
pumped again. Two or more pumping tests are conducted at each well with 
pumping rates which are as different as practicable. For each test, static 


levels are taken and the same procedure is followed in reading piezometers 
as is used in the first test. 


Computation of Transmissibilities 


The layers of the foundation are numbered for convenient reference. The 
numbers increase from the top downward, layer 1 being the top blanket, layer 
2 the top aquifer, etc. The permeability and thickness of a layer are desig- 
nated by kj and 2; respectively, where the subscript indicates the layer. The 
permeabilities considered are those in the vertical direction in blankets and 
those in horizontal directions in aquifers. The vertical transmissibility of a 


k 
blanket is defined to be C; = =, where i must be an odd number. The hori- 


zontal transmissibility of an aquifer is Cj = kjzj , where i must be even. 
Computation of transmissibilities from pumping test results are made by 
the formulas presented below. These formulas are based on the assumption 
that the water table is horizontal and is located above the top blanket. Field 
conditions in the Everglades, where pumping tests for the Central and 
Southern Florida project were made, always satisfied that condition. Some 
modification of the method is required if this condition is not satisfied. For 
example, laboratory tests would be required to obtain the permeability of the 
top blanket where the water level is well below the surface. The contri- 
butions of the various aquifers to the well discharge, Qj where i is even, must 
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be known. Drawdown data should be expressed as drawdown curves—plots of 
drawdown vs. log r, where r is the radial distance from the pump well. In the 
vicinity of the pump well the drawdown curves are straight lines and the total 


radial flow in each aquifer is a constant. From Darcy’s law, the equation for 
radial flow in an aquifer is 


where h, is the drawdown in the aquifer. . 
dhi _ dhj Ologr . dhi_. . jog 
dr dilogn dr ad(logr) 
d(/og r) 
the slope of the drawdown curve. 


Q,(r) = C; Sjcr) 10g 2736; 5,07) 


In the vicinity of the pump well Q; (r) and S; (r) have constant values, Q; and 


Let 


2.73 Si. 


This is the formula used for computing the horizontal transmissibilities of 
aquifers. 

In the case of one aquifer, the vertical transmissibility, Cj, of the blanket 
is determined from the “radius of influence”, re. re is defined as the value 
of r at the point where the extended straight-line portion of the drawdown 


curve intersects the line hg = 0. The ‘radius of influence” satisfies the re- 
lation 


aT 2S 
Since ’ 


the vertical transmissibility, C,, of the blanket can be computed when r, and 
Cy are known. 
In the case of two aquifers, the transmissibilities of the blankets can not be 

computed by a method similar to that used in the case of one aquifer because ’ 
the ‘radii of influence” are not independent of the constants of integration in 

the solution of the differential equations. The vertical transmissibilities of 

the blankets are obtained from the differential equations, the drawdown 

curves, and the horizontal transmissibilities of the aquifers. Two values of : 


r (r, and rg) at which the slopes of the drawdown curves are measurably 
different are selected. 


Qi lr) 2.73 Ci Sir), 


— 
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The C,’s are known and the Sj’s can be measured on the drawdown curves. 
Water can enter the aquifers only by passing through the blankets, so 


AQ, , 


From Darcy’s law, the equations for flow in the blankets are 


Cyr dr. 

rh, dr 


/ 


When values of the integrals are known, these equations can be solved for C,; 
and C3. The values of the integrals may be obtained by numerical integration, 
or by plotting the curves and planimetering the areas. 

The method used for computing transmissibilities of the blankets in the 


two-aquifer case is used in the three-aquifer case also. The formulas re- 
quired are: 


A Q;*2.73 G+246), 
AQ, 
AQ,°4Q,+4Q,, 

AQ,:A4Q, 


Thus, 


A rh, dr , 


r 


AQ, -AQ,- 
| 
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Computation of Underseepage and Uplift Pressures 


The region being considered is a strip of the foundation one unit wide and 
normal to the dam. Qj;, i being an even number, is the rate at which water 
seeps past the center line of the dam in aquifer i. The sign of Q; indicates 
whether the direction of flow is the same as or opposite to the direction arbi- 
trarily chosen as positive on the x axis. Q is the total seepage rate through a 
one unit wide strip of the foundation. Q is always taken to be positive. The 
symbol L is used for the base width of the dam. H is the head differential 
across the dam. The dam is assumed to be impervious, so no water seeps 
through the top blanket in the region beneath the dam. For derivations of all 
equations see Appendices A and B. 

The One-Blanket-One-Aquifer Case.—The underseepage per unit length of 
a dam or levee is given by the formula 

Q: C2 H where 


C / 
ae+f/— ond 
2 


If the soil properties are obtained from a pumping test, then the permeabili- 


ties need not be known since Co and r, are determined directly and “a” can 
be obtained from the relationship 


(r. is defined as the radial distance at which the extended straight-line portion 
of a curve of drawdown versus log r intersects the line of zero drawdown). 
In the portion of the aquifer located beneath the dam the hydraulic head re- 


ferred to tail water is given by 


where X = 0 at the centerline of the dam and X = L/2 at the toe of the dam. 


The Two-Blanket-Two-Aquifer Case.—The following quantities are used in 
the formulas: 


C4 


“2 


) 


— 
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32 
The underseepage per unit length of dam is given by the formula 
Q= -[Q, +Q,/, where 
A Qe YA Q4 
Qa = : 
A 


] 


(3-1) fem 


clo) [meet] 


Z+#l -cfivet]ef g-F 
2(r-i) 


A substantial simplification of the formulas for underseepage without great 
loss of accuracy can be accomplished by making the assumption that no water 
passes through the portion of the lower blanket located beneath the dam. In 
this case 


Q=-[02+9,/ , where 

0, = C2H [#(m-1) , 

ofl Zrgl 
CaH [. #m-l) g(1-J) 


and 
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Since transmissibilities rather than permeabilities are used in the formulas 
for underseepage, the permeabilities need not be known if the transmissibili- 
ties are obtained from a pumping test. The two-aquifer case differs from the 
one-aquifer case in that there are no convenient quantities like x, and re to 
simplify the application of pumping test results to computation of under- 
seepage. 

In the portion of the upper aquifer located beneath the dam, the hydraulic 
ivcad referred to tail water is given by 


Here X = 0 at the centerline of the dam and X = L/2 at the toe of the dam. 
An approximate formula can be obtained by assuming that no water passes 
through the portion of the lower blanket located beneath the dam. 


This formula can be written in another way for comparison with the precise 
formula: 


The approximate formula is conservative in that it gives uplift pressures that 
are too high under the downstream half of the dam. 

The Three-Blanket-Three-Aquifer Case.—Formulas for underseepage that 
take into account the flow through the two lower blankets in the region be- 
neath the dam are too unwieldy for convenient computation since the determi- 
nants involved can not readily be reduced to lower order than the fifth. For 
this reason they are not presented here. It is assumed that no flow through 
the lower blankets occurs in the region beneath the dam. The following 
quantities are used in the formulas: 


C3 Cs _Cs 


Age , 


b = %2 4 %32 (154 
C= %2 


pase are the roots of 2°-az*+bz-c = O, 


P,9,r are the positive sguare roots of pi girs 
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A 5 


Cay 
54 
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The underseepage per unit length of levee is given by the formula 


Q = -[Q2+ Qy + Qe] 


A /2epl "Zegl Zerl 


pA, gaz , 


, where 
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As before, the transmissibilities, rather than the permeabilities, of the layers 
of the foundation are the soil properties required for computation. 

The precise formula for pressure distribution in the top aquifer would be 
cumbersome. The formula presented here is approximate, based on the as- 
sumption that no water flows through the portions of the lower blankets locat- 
ed beneath the dam. In the portion of the top aquifer located beneath the dam, 


the hydraulic head referred to tail water is given approximately by the 
formula 


where X = 0 at the centerline of the dam and X = L/2 at the toe of the dam. 


This formula gives uplift pressures that are too high under the downstream 
half of the dam. 


Uses of Computations 


The computation methods described above were developed for estimating 
the seepage under levees in south Florida where the foundations are composed 
of relatively thick beds of very pervious marine limestone with thin interbeds 
of dense fresh-water limestone. The methods are also applicable to seepage 
under any other water-retaining structure on a foundation consisting of alter- 
nating impervious and pervious materials. 
Underseepage computations are particularly valuable in cases where water 
conservation, either for water supply or for the generation of power, is a 
primary consideration. Such computations can indicate whether or not the 
cost of the provision of an upstream blanket or a cut-off can be justified by 
the benefits which will be realized. All too often such features are incorpo- 
rated in designs without a complete analysis of their economic justification. 
Computations of uplift pressures are very useful in indicating whether or not 
pressure-relief devices such as toe drains, drainage blankets, and relief wells 
are required to prevent possible piping. 
It is not meant to imply that exact, or even very accurate, predictions of 
underseepage quantities and uplift pressures can be made by the methods de- 
scribed above. The methods require an idealizing of the non-uniform con- 
ditions actually encountered in nature. Furthermore, it is unlikely that the . 
basic data used in the computations can be obtained with a high degree of ac- 
curacy. It should be noted, however, that these qualifications are equally 
applicable to all soil mechanics design computations. 


CONCLUSIONS 


' Methods for computing seepage and head distribution through foundations 
made up of alternating impervious and pervious layers have been developed. 
These methods can be used for determining seepage quantities and uplift 


pressures in systems which are too complex for the practical application of 
flow-net analyses. 
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SEEPAGE PREDICTION 
APPENDIX A 
Derivation and Solution of Differential Equations for Head Distribution 
in Flow Through Layered Systems 
Derivations of the Differential Equations 


The One-Aquifer Case 


The differential equation for flow in a single blanketed aquifer can be de- 
rived by considering a region of the aquifer with length dx, width dy and height 
Z9. 


The total flow out through the sides of the region is equal to the total flow 
in through the top. 


where V, and V, are velocity components, p is the hydraulic head in the 
aquifer, andp, is the constant hydraulic head at the top of the blanket. The 
last term in the above equation was obtained from Darcy’s law, V = -kVp. 
Dividing the equation by dx . dy gives 


_k 


Let h=p.-p, C, Darcys law 


Op 
Substituting in the equation, 
Oh 


| 
oe. 
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an, 
“Oy? 


This is the required differential equation. In cylindrical coordinates the 
differential equation is 
r $2 


Multiplying by r2 gives 


The Two-Aquifer Case 


The differential equations for the two-aquifer case are derived by the same 
method as was used for the one-aquifer case. Consider the regions of the two 
aquifers inclosed by a prism with length dx and width dy. For each region 


the total flow out through the sides is equal to the total flow in through the top 
and bottom. 


dx z,dy+ dy z2dz “ky dzxdy 
/ 
dzdy, 
AX zy dy + Wile cy ay dz +k, dzdy 


From Darcy’s 


Let 


7 
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Substituting in the equations, 
2 
be Che- C3 (hy-he) J 
0% 


Let 


Ci Ca C3 


12 Ce 32 C4 


2 


This is the required pair of simultaneous differential equations. In cylindrical 
coordinates 


or 
Multiplying by r2 gives 
Oth; dhe , d*hg 


2 2 


The Three-Aquifer Case 


Derivation of the differential equations for the three aquifer case is so 
similar to that for the two-aquifer case that there is no need for including the 
details of the derivation. The simultaneous differential equations are 


2 2 
od*ha 


ox? dy? = O54 (ha-he)- %5y 


4 
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z 2 


Solutions of the Differential Equations 


The One-Aquifer Case 


Sheet Flow. —Consider the flow toward a dam in a strip of the foundation 
one unit wide. Let x equal zero at the heel of the dam and increase in the up- - 
stream direction. Only the region x >0 is considered here. Since h = h (x), 


By2 = 0 and the differential equation becomes 

. Pe 

h=0. 

(D-a)(D +a) 
The general solution of the differential equation is 
h-=M,e% , 


Where M, and Mg are arbitrary constants. From the physical conditions of 
the problem, 


h=0. 


co 


Consider the tangent to the curve h versus x at x = 0. The x-coordinate, x,, 
of the point where the tangent intersects the line h = 0, is called the “re- 
sistance” of the blanketed aquifer. 


| 
Since q 
a as M,=0. 
h(t) =-aM, 
=-@ M, --ah(o). 
| 
/ 
Ep = 
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Radial Flow. —Consider the flow toward a well with its centerline at r = 0. 
2 
Since h = h (r), a = 0 and the differential equation becomes 
d? A, dh gtr? A. 
dré 


6-rD. 
r? +52 (6-5 0%. 


-@ t,2) h = 


I, (ar) and Kp (ar)—the zero-order modified Bessel functions of the first 
kind and second kind respectively—are linearly independent particular so- 
lutions of this differential equation. The infinite series expressions for these 


functions can be found in Tables of Integrals and Other Mathematical Data, by 
Dwight (formulas 813.1 and 815.1). They are: 
L 


Zn 


where Y is Euler’s constant, 0.5772... . Values of the function 


; Hin) = K, 


are tabulated on pages 236-242 of Tables of Functions, by Jahnke and Emde. 
The general solution of the differential equation 


-@ *,*) hy = O 


is h = NylIp (ar) + NgKy (ar), where Ny and Ng are arbitrary constants. From 
the physical conditions of the problem, 


lim h= 


co 


as and tim K,(ar)-0, 
h = N, K, (ar) 


When r is small, then, approximately, 


h=-N,(7+Iln 


| 

| 
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The approximation gets better as r becomes smaller. The drawdown 
curves (plots of h vs. log r) obtained from observations made during pumping 
tests exhibit straight-line portions for small values of r. 
re, the value of r at the point where the extended straight-line portion of a 
drawdown curve intersects the line h = 0, is called the “radius of influence” 
of the blanketed aquifer. ‘ 


-Nz (7+#/n 0. 


The Two-Aquifer Case 


Sheet Flow.—Consider the flow toward a dam in a strip of the foundation 
one unit wide. Let x equal zero at the heel of the dam and increase in the up- 
stream direction. Only the region x 20 is considered here. Since hg and hq 


2 
are independent of y, 352 - = 0 and the differential equations become 


(h4-he) = 0. 


[D?- 12+ %2)/ he 


2 
+/D sy/hy+0. 
This pair of simultaneous equations will be solved in the following manner. 
The variable hq will be eliminated between these equations. The resulting 
equation will be solved for hg. When hg is known, hq will be obtained from 
the first equation of the pair. Operate on the first equation with [2 . ot34] 
and then subtract the product of a@39 and the second equation. 


LD*-*s4] [D* #%30)]h, -%32%u hr, - 
[D* ("12+ sy) D*+ h, 


> 


[D4 aD*+b] h, 


Let 
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Let 


he=0- 


The general solution of this differential equation is 
x 


where the M’s are arbitrary constants. Now that ho is known hq can be ob- 
tained from the equation 


2 
o 
32 32 
ha Mz et +mM,e4 +mMye9* , 
The general solution of the pair of differential equations is 
h,=Me + M,e7*+M,e™ , 
9% 


where the M’s are arbitrary constants. 
From the physical conditions of the problem, 


lim he /im hy = 0 


= 00 

et ic and as 

>» 


h,: /M,et@+m 


Since 
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The quantities corresponding to the “resistance of the blanketed aquifer in 
the one-aquifer case will now be obtained. 


h(0)=-(f Mp + 


-(4M, +9 My) =- 


For the upper aquifer, 


Mz+Ma 
Similarly, for the lower aquifer, 
Mzt+mMa 
M4 
In the region beneath an impermeable dam no water can pass through the 
top blanket. The derivation of differential equations already presented ap- 


plies to this case when the term associated with flow through the top blanket 
is omitted. The required differential equations are 


= 


dy* 
/h, + hg =O 
h, +[D* /h, - 0 


Solution for hg can be effected by eliminating hg between these equations. 


he first equation of the pair then gives h4. Operate on the first equation with 
D2 - asa] and then subtract the product of og and the second equation. 


54 *52/h, -%2 “se h,=0. 


Er 

=; 

Since 

| 
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[D*- (32 + 34) D*] h, 


« 


D*[D-c][ D+c] hy ‘ 


h, , 


where the R’s are arbitrary constants, 


S 


The general solution of the pair of differential equations is 


cz 
27k, +kse +Rae 
Cx “Cz 


where the R’s are arbitrary constants. 
Radial Flow.—Consider the flow toward a well with its centerline at r = 0. 


2 
Othe oat ait = 0, and the differential equations become 


%32) r ‘lh, +%z 
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d=rD. 


[8* (ie hg =O 


Analogy with the one-aquifer case indicates that the general solution of this 
Be of equations should be 


(fr) +N, K, (fr) +N, L, GDN, K, » 


where the N’s are arbitrary constants. This solution can be verified by 
substitution in the differential equations. The required formulas for deriva- 


tives of Bessel functions, which may be found in Tables of Integrals and Other 
Mathematical Data by Dwight, are 


“ak, =-nK,-xK,_,,and K,=- 
From the physical conditions of the problem 


/im h, =/im fig =O. 


—oo 


Since 


and lim K,(%)- 0, 
N, 
h,=N,K, (tr) + NgK, (ar) , 
hg> N,Ky (fr) +r Ng &, (9r)- 


The drawdown curves obtained from observations made during pumping 
tests have stright-line portions near the pump well. When r is small, then, 


approximately 
K, (fr) 


The equations of the straight-line portions of the drawdown curves are 


| 
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hy=-({N,tm N, In Inf - 
The quantities corresponding to the “radius of influence” in the one-aquifer 
case are the values of r at the points of intersection of those lines with the 


lines hg = 0 and hq = 0, respectively. 


The Three-Aquifer Case 


Sheet Flow.—Consider the flow toward a dam in a strip of the foundation 


+ %o)/hy hg =O, 
34/2 +[D*- +54 Ne =O) 


< 5 | Ne 0. 
This system of simultaneous equations will be solved in the following manner. 
hg and hq will be eliminated to give a differential equation that can be solved 
for hg. Then hq will be obtained from the first of the ‘[o a and hg 


from the second. Operate on the second equation with | D2 - Os | and sub- 
tract the product of a54 and the third equation. 


Multiply this eqyation by -g9, then operate on the first of the differential 
equations with |D4 - (cgq4 + 054 + 56) D2 + O34 56 | and add results. 


x 
| 
one unit wide. Let x represent distance from the heel of the dam. Since hg, fe 
4 and hg are independent of y, an e eren 
equations are 
i 
| 
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%32 “x, /h, = O. 


(%5y 56) + “se D2 fi,:0. 


a= % 2 +%54 
b = %2(%syt%5u 56) + (%54 50) 


c= “su “Se 
[D°-aD*+bD +0. 
Since the physical conditions of the problem and the solutions in the cases of 


one and two aquifers indicate that there is nothing periodic about hg, the above 
operator should factor as follows: 


[0° aD + bD* 


where p2, q2 and r2, the roots of #3 - az? + bz - c = 0, are real and positive. 
Pp, q, and r are taken to be positive. 


h, 


The general solution of this differential equation is 
hi, =M, +M, e°"4M,eF* +M,e *+M,e 


where the M’s are arbitrary constants. Now that hg is known, hq and hg can 
be obtained from the first two of the simultaneous equations. 


hy = [ig + %52)-D*] hy, 


Let 


Z 
%32 


Let 
| 
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Let 

54 

The general solution of the system of simultaneous differential equations is 


where the M’s are arbitrary constants. From the physical conditions of the 
problem, 


lim h,=lim hy lim h, =O. 


rz 


h, = + , 


e 


+nM,e"* . 
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APPENDIX B 


Derivations of Formulas for Computation of Seepage and Uplift 
Pressures in Layered Systems 


One-Aquifer Case 


Seepage Under a Dam 


The significance of the “resistance”, x,, of a blanketed aquifer is that, as . 
far as h (0) and Q (0) are concerned, the infinite-blanketed aquifer in the 
region x > 0 is equivalent to an aquifer with the same permeability and thick- 
ness extending from x = 0 to x = xy and so arranged that the upper and lower 
boundaries are completely impervious and h(x,) = 0. If there is no break in 
the blanket, then, for the purposes of computation, the infinite-blanketed 
aquifer on the pool side of the dam may be replaced by an aquifer of length x, 
as indicated above. By symmetry, the same thing may be done on the other 
side of the dam. Let L be the base width of the dam and H be the head differ - 
ential. The underseepage, Q, per unit length of dam is given by the formula 


Ole 


where 


Uplift Pressures 


In the portion of the aquifer beneath a dam the hydraulic gradient is a 
constant equal to Q/Co, so the hydraulic head varies linearly with x. Let x 
equal zero at the center line of the dam and increase in the upstream di- 
rection. The curve h vs. x should be symmetrical with respect to the point 


(O,H/2). 
hea Gy, +5) ; 
Zz 2 2 
where h is the absolute value of the hydraulic head referred to head water. 


By symmetry, the same equation applies when the direction of x is reversed , 
and h represents hydraulic head referred to tail water. 


Two-Aquifer Case 


Seepage Under a Dam 


In the two-aquifer case, the quantities corresponding to the “resistance” of 
a blanketed aquifer have values that depend on boundary conditions as well as 
the constants of the formation. For this reason it is not convenient to use 
them in the derivation of formulas for underseepage. Let H be the head differ- 
ential, and let x = 0 at the heel of the dam, x = -L/2 at the centerline of 


q 
| 
I | 
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the dam. Each of the curves (hg vs. x and hq vs. x) should be symmetrical 


with respect to the point (-L/2, H/2). The underseepage, Q, is the flow under 


the dam, 


Expressions have already been obtained for hg(x) and h4(x) in the regions 


When 


-lp<x<0, 
ho =R, + , 
hy +Rox%+ sRye~* 


When 


hyo = +M, 


Considerations of symmetry and of continuity of head and gradient give a set 


of simujtaneous equations that can be solved for Qo and Q4: 


he (-Ye)=hy(-Y2)= Yo 
lint = lim, (6-24). 


From equation (a), . 


(a) 


(b) 
(c) | 
(a) 
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From equation (b), 


R, +cRze 


Cc 
+scR3e %-scRye” - 
From equation (c), 
R, +R #Ry = =-Mz +My 
R, = = + +mMy 


From equation (d), 


R, =-*fMz-9My » 
Rot scks -s¢Ry=-JtM, 
Rewriting the above equations, 
cL -0 
2 kag +e2 Ry 
R,-5 Re 
cl cl Q 
=O, 
th 
R,+sce 2 R3-sce * hy 
Ry -M, -M, =O, 
A, SRy-yMz-mMy =O, 
R, + - +#M,+9M, 
R, + sch; -scRy +jfM, +mgMy = 
Cramer’s rule will be used to solve these equations for Qo and Q4. 
et 
Z 
/ -+= O O 4 O 
C CL 
scé2 -sce? 0 -fy\s 


cl 
-cRye * ? 
| 
| 
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S -/7) 
-C f 
/ SC -sc sf mg 


o(j- s)(1+e° gfm-s) 


/ 
0 
/ 
/ 
O 


c(e% 


/ 


? 
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/ C..@ 
A- 
which reduces to 
: 
C2 Cy 
| 
Zz -Cé 
Ko scé 2 -sce? 0 0 - Cy 
which reduces to 
_cl cl 
o 6-ce2 O 0 
cl 
lo 0 0 O 
“J 
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/ 


-m 0 O 


Jt# mg 0 o 


which reduces to 
c(s-1)(14 
The underseepage, Q, per unit length of dam is given by the formula 


Q=-/@, #Q, /, where 
Qo - 


-Ng 


cl 
CW | 6-Ni-e%) 


9(/-m) 


AQy CyH 


/-m 


/ SC -SC 
| 
g 
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Li-i) 

/-/ g(n-j) 


The above expressions for AQ9, AQ4, and A were multiplied by 2 CoC4 e a 


Since only the ratios of those quantities are used, the results are not affected. 
Approximate formulas for computing underseepage in the two-aquifer case 
can be obtained by assuming that there is no vertical flow through either 
blanket in the region beneath the dam. 
When 


ho = M, etx My e 


hy +mMyé 


-§<x<0; 


? (i= 2,4) 


dh 
since Qj = Cj = = a constant, and hy (-L/2) = H/2. Continuity of head and 
gradient in each aquifer requires that 


hitk) = Jim, » (i= 2,4), 
lim. Dhi(x) = fim Dh, (x), ((-2,4). « 


From equation (e), 


(e) 


ASCE 1 | 
| 
When 
H LQz = M. 
49 ce Mp 
H LQ¢_/;M My, 
2C4 
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From equation (f), 
Qz = - Mz - 9 
2 

(g) 
=~ -gmMy - 


Multiplying the last pair of equations by L/2 and subtracting from the pre- 
ceding pair of equations gives 


= Mz + (1495) Mo > 
- tm (1+ 95) Mz 


Ali ml? 4 


By Cramer’s rule, 


Substituting these values of My and Mg in equations (g) gives 


2+fL * 
H [- #(m-/) 
IF 


The underseepage, Q, per unit length dam is given by the formula 


where Qo and Q, are obtained from the preceding pair of equations, 


Uplift Pressures 


In the portion of the upper aquifer beneath a dam the absolute value of the 
hydraulic head referred to head water is given by 


: 
| 
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cz 
R, + Ree O), 
where x = 0 at the heel of the dam. The first four of the eight simultaneous 
equations that were used in determining underseepage can be solved for the 
coefficients in this equation. 
cl 
e= Ry 


cl Q 
- Fey Cy 
cl 
O 2 e 2 
/ -5 + 
A = cl y 
/ -Ce a 
cl 
/ -sce 2 
which reduces to 
A=2c(/-s)-> 
Cz 
Qa / cl 
Cy sce -~sce 
L | Ge_ ]. 
cH (/-s)4e s / 
Since 
s=- £2 and Q=-Qe- 
Cy 
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H 
2 

Q 


2 


> 
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AR, = cH Cy 
cL 
0 @ * 
AR, = - SE 
Cz cl 
0 Q4 sce 2 -sce 2 
C4 
cl 
H 
BR; = cl 
/ AZ -ce 
C2 
Qe 
/ Ce 
( gz ge) 
C2 Cy 
/ 0 
0 / ce # H . 
a 
0 / sce = = | 
(4 2 
C2 Cy 
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Using Cramer’s rule, 


Xe cX , 
[02 Cy/ 


Since the curve hy versus X should be symmetrical with respect to the point 
(O, H/2), this equation also applies in the interval ( + <X< 0). Thus, 


Qe_Qe/sinh cX (Leyek 


Because of symmetry this equation is also correct when X increases in the 
downstream direction and hg is the hydraulic head referred to tail water. 

An approximate formula for the head in the portion of the upper aquifer 
that is beneath the dam can be obtained by assuming that the part of the lower 


- 

2= Cy | | 
ol 
2c(/-s) / C2 _Cy | 
2°2 Untly) CotCy 20 
& | 
Cotly (2 ia) 
Let 
| 
_H L i 
| 
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blanket that is beneath the dam is impervious. In this case the formula de- 
rived for a single aquifer will apply when h and Q are replaced by hg and -Qo, 
respectively. 


($sxs$)- 


This formula can be written in another form for comparison with the precise 


formula. 
, Wet Qy Q2tQy 


Bethe CaQe-C2 Qu . C2 Cy 
Ce(Cet# Cy) C2 Cy 


Qy/_ Ce 
Co Cy/ Cotly 


Q2tQy - Cy 
C2 tly -Cyst ly 


#Cy 


Cy /-S$ 


L 
Tree 


Three-Aquifer Case 


Seepage Under a Dam 


The first method presented for computation of underseepage in the case of 
two aquifers is too cumbersome to be practicable in the three-aquifer case. 
Approximate formulas for computing underseepage in this case can be ob- 
tained by assuming that there is no vertical flow through any blanket in the 
region beneath the dam. Let x = 0 at the heel of the dam, x = -L/2 at the 
center line of the dam. Each of the curves h, versus x (i = 2, 4, 6) should be 
symmetrical with respect to the point (-L/2, ‘H/2), where H is the head 
differential across the dam. 


| 
| 
+Qy=-Q 
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+ Mee” » 
hy=Mp + gMye + > 
h, =kMpe” +mMye + 


dhj 
since Qi = Cj Gy = a constant, and hj (->) = 3: Continuity of head and 
gradient in each aquifer requires that 


lim. hj la)= lim, 
Le O 


lim (x) = lim, Dhi (6 =2,4,6)- 


From equation (h), 


M. + Ma 


Re = k Mo 

From equation (i), 
Qe 
Cz 
at =-fpMz -99My >» 


C, =~kpM,z -mQMy, -nrM, 


=-pM, - gMy -rM, 


ASCE 
When 
When 
: -5exs O ; 

hi = +3 » (i 16) > 
f+ 
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Multiplying the last three equations by L/2 and subtracting from the preceding 
three equations gives 


(1405 +(/+9 5) My Me 


By Cramer’s rule, 


(110 $4 


HA, __ Hz 


Substituting these values of Mg, Mg, and Mg in equations (j) gives 


* 


*92+gL */2FrL 


CoH k pA, 4m rls 
A OL 27QL 
The underseepage, Q, per unit length of dam is given by the formula 


Q = - [Ay + +Q, / ? 


where Qo, Q4, and Q¢ are obtained from the preceding three equations. 


Let 

A=|f 9 J|, A,=\/ 9 / \,hg= 

kmn / mn k/n 

\ 
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Uplift Pressures 


The first method presented for computation of uplift pressures under a 
dam in the case of two aquifers would be very cumbersome in the three- 
aquifer case. An approximate formula can be obtained by assuming that there 
is no vertical flow through the part of the middle blanket that lies beneath the 
dam. In this case the formula derived for a single aquifer will apply when h 
and Q are replaced by hg and -Qg, respectively. 
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DESIGN AND PERFORMANCE OF VERMILION DAM, CALIFORNIA 


K. Terzaghi,! Hon. M. ASCE and T. M. Leps,2 M. ASCE 
(Proc. Paper 1728) 


SYNOPSIS 


This paper describes the design and performance of a zoned earth dam 
located in the Sierra Nevada, California, on marginal glacial deposits. The 
pattern of stratification is so intricate that it was impracticable to determine 
all the essential geological features of the subsoil in advance of construction. 
Therefore, it was necessary to check on the design assumption during con- 
struction and to modify some details of the original design in accordance with 
the findings. As a result of this procedure the finished structure is as safe 
and satisfactory as if all the essential properties of the site had been known in 
advance of construction. 


INTRODUCTION 


The damsite is located at an elevation of about 7500 feet above sea level, 
in the valley of Mono Creek, a northern tributary of the South Fork of the San 
Joaquin River in central California. The site was investigated for the first 
time in 1918 to determine the watershed and capacity of the reservoir. Ex- 
ploratory drilling in 1924 revealed unfavorable foundation conditions for the 
concrete dam then being considered, and the project was temporarily aban- 
doned. It was reconsidered in 1940 on the basis that an earthfill dam could be 
shown to be feasible. The earthfill dam, as built, has a maximum height of 
160 feet above the original ground surface, and the crest has a length of about 
4,300 feet. Construction was started in May 1953 and the dam was completed 
at the end of 1954. The reservoir formed by the dam adds 125,000 acre feet 
to the total storage capacity of the reservoirs which have previously been 
built by the Southern California Edison Company in the San Joaquin basin. A 
Note: Discussion open until January 1, 1959. To extend the closing date one month, 

a written request must be filed with the Executive Secretary, ASCE. Paper 1728 is 

part of the copyrighted Journal of the Soil Mechanics and Foundations Division, 


Proceedings of the American Society of Civil Engineers, Vol. 84, No. SM 3, 
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description of the entire hydroelectric power development utilizing the water 
resources of the basin, and a brief account of the construction of the 
Vermilion Dam have been presented by Mr. R. W. Spencer, Manager of Engi- 
neering Department, Southern California Edison Company, in the March 1955 
issue of Civil Engineering. 


Geology of the Damsite 


Before the Pleistocene ice age, the present location of Mono Creek Valley 
was occupied by an erosion valley carved out of granite. Locally the surface 
of the granite was covered by the remnants of a flow of pre-Pleistocene 
basaltic lava. At the damsite, the bottom of the old riverbed is located at a 
depth of about 250 feet below the deepest point of the present valley. An out- 
crop of the basaltic lava is located a short distance upstream from the dam- 
site, west of the old course of Mono Creek. 

During the Pleistocene ice age, the valleys descending from the summit re- 
gion of the Sierra Nevada were repeatedly invaded by valley glaciers descend- 
ing to an altitude of a few hundred feet below that of the damsite. The glaciers 
deposited subglacial, marginal and fluvio-glacial materials in the valleys. At 
the damsite, the youngest member of the glacial deposits consists of a set of 
low, terminal moraines resting on a locally very impervious glacial till. At 
the west end of the damsite, the remnants of a kame terrace were encountered. 

The upper portion of the fluvio-glacial sediments located beneath the young- 
est till sheet contains one large and several smaller lenses of silt with a very 
variable thickness. The continuity of at least the youngest fluvio-glacial and 
glacial lake deposits was disrupted by small erosion valleys carved out by 
streamlets emerging from the glaciers and subsequently backfilled with silty 
to clean sand and sand-gravel mixtures. Therefore, the pattern of stratifica- 
tion of the subsoil at the damsite is so erratic that it was impracticable to 
determine its average coefficient of permeability on the basis of the results 
of laboratory tests. 

Upstream from the damsite, the subsoil consists of coarse-grained, some- 
what silty, fluvio-glacial and post-glacial river deposits with an almost hori- 
zontal surface, underlying an area of the valley floor with a length of several 
miles and a width of several thousand feet. Formerly, Mono Creek meandered 
on this floor. In recent geological times, the meandering creek started to 
erode, forming incised meanders with steep side slopes. At the damsite, the 
creek has cut to a depth of about 40 feet below the old valley floor. 

One of the important considerations indicating the reasonableness of con- 
structing a dam on such a site was that the Sierra Nevada Mountains contain 
nearby many natural lakes and tarns formed entirely by glacial debris fills 
left across their outlets. These have stood for centuries without important 
water loss from percolation. Why then could not an artificial fill be made 
to approximate these natural fills and thereby obtain a sound structure, built 
entirely of country materials? 

The details of the glacial geology of the damsite were investigated by 
Mr. J. H. Birman during the summers of 1950 and 1951, under the general 
supervision of Dr. Ian Campbell of California Institute of Technology. The lo- 
cation of the terminal moraines and the kame terrace at the west abutment 
are shown in Figure 1. Figure 2 is an idealized geological profile of the sub- 
soil along the centerline of the dam. The results of the geological survey 
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furnished the basis for a tentative interpretation of the boring records and for 
making reasonably realistic assumptions concerning the subsoil conditions 
for design purposes. 


Subsoil Explorations 


In 1924, 18 boreholes were made at and in the proximity of the present site 
of the dam. On the basis of the results, the project was temporarily aban- 
doned. In 1946, the site exploration was resumed under the general direction 
of R. W. Spencer and under the direct supervision of the co-author, T. M. 
Leps. The exploratory work included 14 borings with an average depth of 100 
feet and a maximum depth of 269 feet. The exploration also included the dig- 
ging of 2 test pits with a depth of about 30 feet each, of 3 open cuts across 
terminal moraines and several bulldozer trenches on the steep slopes of the 
present erosion valley of Mono Creek. Eight of the 14 drill holes were located 
along the centerline of the dam and spaced abouy 600 feet apart. The other 
holes were drilled in the proximity of the toes of the dam. 

The crest of the dam is about 4300 feet long and the base has a maximum 
width of 800 feet. If the drill holes had been spaced about 20 feet both ways, 
subsoil exploration would have required several thousand drill holes and a 
prohibitive amount of testing. Yet, judging from the pattern of stratification 
visible on the sides of open cuts, for instance on the exposure shown in 
Figure 3, the results would still have left a wide margin for interpretation. 
Therefore, at a site like that of Vermilion Dam, the exploratory drill holes 
merely serve the purpose of supplementing the results of the geological sur- 
vey. At the site of Vermilion Dam, 14 borings sufficied to obtain all of the 
required information. Any additional exploratory work would have increased 
the cost of exploration without furnishing any significant contribution to our 
knowledge of the site. 

The tests on the samples recovered from drill holes and test pits included 
grain size determination and permeability tests. Even within the fluvio- 
glacial deposits the value of the coefficient of permeability k was extremely 
variable. It ranged from between 10-1 and 10-6 cm/sec. On account of the 
silt content of the sediments, the average was rather low, and the coefficient 
of permeability in a vertical direction, ky, was consistently smaller than the 
corresponding value kh, for horizontal directions. 

The results of the permeability tests made it possible to obtain at least an 
upper and lower limiting value for the loss of water, Q, due to seepage out of 
the reservoir, using the flow net method. Assuming: 


1. ky, = 9 x 10-3 cm/sec and ky = 10-3 cm/sec, Q = 29 cfs and 
2. ky, = ky x 10-3 cm/sec, Q = 4 cfs 


The observations described under a later heading entitled “Seepage Losses” 
will show that the total leakage at full reservoir is roughly equal to 6 cfs 
which is slightly greater than the foregoing estimated lower limiting value. 


Construction Materials 


Within the area surrounding the damsite, the following construction mate- 
rials were available in significant quantities: stony till, fluvio-glacial sand 
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and gravel, silty pond deposits, and the weathered top layer covering the un- 
weathered river deposits upstream from the damsite. Clay deposits were 
conspicuously absent. Even the pond or lake deposits contained no more than 
a trace of clay. This is probably due to the steep gradient of the water 
courses which were fed by the glaciers. Most of the clay was removed in 
suspension. 

Because of the absence of clay, all efforts were concentrated, in the earlier 
stages of the exploratory work, on locating suitable materials for the ‘im- 
pervious” section of the proposed earth dam and for the construction of an 
“impervious” blanket to be placed upstream from the dam. After the borrow 
pits were opened, it was found to be equally difficult to find very pervious ma- 
terial, because the bulk of the glacial outwash contained at least traces of silt 
or rock flour. 

On the basis of the results of the geological survey, four general areas 
were selected as potential sources for suitable construction materials. The 
location of these areas with reference to the damsite are shown in Figure 4. 
Area #1, at a distance of about one mile west of the dam consisted of a 
marshy meadow covering the surface of a stratified to varved deposit of fine 
Silty sand and silt which was laid down in a pond. In order to estimate the 
available quantity, auger borings with an average depth of 20 feet and a maxi- 
mum depth of 90 feet, spaced about 400' both ways were made. It was found 
that the laboratory coefficient of permeability of the material which ranged 
between 10-6 and 10-2 cm/sec, averaged 10-5 cm/sec. The natural water 
content of the material was above optimum. Therefore, it was decided to 
lower the water table beneath the meadow about 5 feet by means of drainage 
trenches with a depth of about 10 feet (see Figure 4). Since the deposit con- 
tained more or less continuous layers of fine, clean sand, drainage was effec- 
tive. Most of the material derived from Area #1 was classified as impervious 
and was used for constructing the blanket and core. 

Borrow Areas #2 and #3 were close to and north of the eastern portion of 
the damsite. The water table was located at an average depth of 20 feet below 
the original ground surface. Most of the material encountered in these areas 
consisted of fluvio-glacial sand and gravel, containing scattered boulders with 
a diameter up to several feet. The coefficient of permeability, k, of the mate- 
rial ranged between 10-2 and 10-5 cm/sec, depending on the silt content. The 
most pervious material was placed in contact with the downstream face of the 
“impervious” core and on the base of the downstream portion of the dam. The 
remainder of the downstream portion and the upstream portion of the dam 
were made of the siltier varieties of the borrow material, with a permeability 
range of 10-3 to 10-4 cm/sec. 

Area #4 was located on the wide and almost horizontal floor of the Mono 
Creek Valley upstream and within a distance of from two to four miles from 
the damsite. The deposit was explored by means of auger holes, 10 feet deep, 
spaced 700 feet in the direction of the valley and 500 feet at right angles to it. 
The borings showed that the area was underlain by pervious and semi-pervious 
fluvio-glacial deposits and river sand and gravel. The uppermost 3 ft. of the 
deposit were slightly weathered and relatively impervious. The water table 
was located at an average depth of 9 feet below the valley floor. By excavat- 
ing the uppermost portion of the deposit to an average depth of 5 feet and mix- 
ing the pervious and impervious materials present within this layer, 1.3 mil- 
lion cu. yds. of compacted fill material with an average permeability of 10-5 
cm/sec obtained. Such a blend of the top 5 feet of Area #4 furnished the bulk 
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of the “impervious” core of the dam. The mixing was satisfactorily per- 
formed by excavating the material along a vertical face with a height of about 
5 feet, using Sierra Loaders. A portable sprinkler system was used to add 
moisture to the surface materials. The output averaged 20 cu. yd. per minute 
for each loader. 

The material iocated beneath the 5 foot thick excavated layer was relative- 
ly pervious sand and gravel having an average coefficient of permeability of 
10-3 cm/sec. Some 1.5 million cu. yds. of it was excavated and built into the 
upstream and downstream portions of the dam. Some of the borrow pit exca- 
vation for such material extended well below groundwater table, necessitating 
the excavation of the drainage ditch shown on Figure 4. The pervious materi- 
als were loaded by DW-20 scrapers assisted by two push-cats. 


Design of the Dam 


The design of the dam was governed by the following considerations. 

(1) Due to the scarcity of “impervious” materials, the construction of a zoned 
earth dam was indicated. (2) The average permeability of the substrata is 
high, and (3) the cost of a cutoff of over 200 feet to bedrock would be prohibi- 
tive. Therefore, it was necessary (1) to cover the ground surface upstream 
from the dam with an impervious blanket, and (2) to provide adequate drainage 
facilities along the downstream toe to eliminate the possibility of a failure due 
to piping. 

Figure 5 is a plan of the dam, Figure 6a is a vertical section through the 
highest portion of the dam at about Station 21+00 and Figure 6b a section at 
about Station 28+00. The “impervious” core of the dam rises above the down- 
stream fringe of the “impervious” blanket which covers the natural ground to 
some distance upstream from the dam. According to the original design, the 
fill upstream from the core was to be made of silty sand and gravel with an 
average coefficient of permeability k of less than 2 x 10-3 cm/sec, and the 
fill downstream from the core out of relatively clean sand and gravel with a k 
value of more than 2 x 10-3 cm/sec. 

The water entering the pervious downstream section of the dam through its 
base flows into a toe drain installed in a trench with a maximum depth of 20 
feet. Figure 7 is a typical section through the toe drain. Between Station 22 
and 28 where the natural ground rises beneath the downstream slope to an 
elevation of about 70 feet above the toe trench (see Figure 6b) the toe drain 
was supplemented by an interior drain located between the core and the upper 
edge of the buried north slope of Mono Creek (see also Figure 5). 

Because of the predominantly horizontal stratification of the subsoil and 
the presence in the subsoil of lenses and layers of silt with unknown horizon- 
tal dimensions underlain by very pervious sediments, it was conceivable that 
the hydrostatic pressure acting on the base of silty layers might lift the floor 
of the valley downstream from the dam. Therefore, it was decided to supple- 
ment the toe drain with relief wells as shown in Figures 7 and 10. At the 
maximum section of the dam, where the danger of heave was greatest, the 
wells were drilled to an average depth of 70 feet and the spacing ranged be- 
tween 50 feet and 100 feet. Beyond the maximum section of the dam and the 
outer boundaries of the Mono Creek stream bed, the wells were drilled to a 
depth of about 50 feet and the spacing was increased to 200 feet (see Figure 10). 

The stability computation for the dam was made by the Swedish Circle 
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method. In any investigated section, a large portion of the critical circle was 
located in the natural ground, beneath the base of the fill. The factor of 
safety with respect to sliding along the potential surface of sliding depends to 
a large extent on the average intensity of the porewater pressures along the 
surface. Because of the relatively high permeability and low compressibility 
of the core material, it was justified to assume that the construction porewa- 
ter pressures would be negligible. The porewater pressures due to leakage 
out of the reservoir are determined by the pattern of seepage through the sub- 
soil. Such information was unknown until the reservoir was filled for the 
first time. Therefore, the stability computations were made on the basis of 
what appeared to be the highest porewater pressures compatible with general 
geological characteristics of the subsoil. 

Trial computations showed that the loss of water due to seepage out of the 
reservoir could turn out to be prohibitive unless the core of the dam was sup- 
plemented by an “impervious” upstream blanket. Due to the presence of 
lenses and layers of silt in the subsoil, the foundation of the dam represents 
an irregularly stratified, layered system. Mono Creek and its western tribu- 
tary, Cold Creek, have cut trenches into the system to a depth of about 40 feet, 
thus permitting the water to enter the highly pervious strata located between 
silty ones. Therefore, it was decided during the preliminary investigation to 
seal the bed and the slopes of Mono Creek by means of a blanket extending to 
a distance of about 1300 feet upstream from the core of the dam, and to seal 
in a similar manner the bottom and the south slope of Cold Creek between 
Mono Creek and the west abutment of the dam. 

On the eastern slope of Mono Creek, a silt stratum with a variable thick- 
ness was exposed. The outcrop appeared to extend from a point about 1000 
feet downstream from the damsite to a distance of more than 2000 feet up- 

, Stream from it. The surface of the silt stratum was located at a depth rang- 

ing from 10 to 25 feet. When the preliminary design was prepared, it was be- 
lieved that the silt stratum underlies the entire eastern half of the damsite. 
Therefore, it was proposed to construct east of Mono Creek, a cutoff extend- 
ing from the base of the core into the silt stratum as shown in Figure 8a. 
Subsequently, the modification shown in Figure 8b was proposed to offset ex- 
pected seepage through the east abutment area. However, before the con- 
struction drawings were prepared, exploration disclosed that the dam was 
located above a narrow but deep pervious gap in the silt stratum. Therefore, 
the location for the cutoff was shifted into the position shown in Figure 8c. 
Core and cutoff are connected with each other by an “impervous blanket.” 
The thickness of the blanket is 13 feet from the core to the cutoff, whereas 
the thickness of the blanket in the valley of Mono Creek decreases from 
about 13 feet at the core to 6 feet at its upstream boundary. 

The service spillway with a discharge capacity of 1200 cfs is located at the 
east abutment of the dam (see Figure 5). It is a concrete weir with a crest 
elevation of 7634.5, equipped with one 15 ft. x 8 ft. radial gate. Floods in ex- 
cess of 1200 cfs can escape over a broad emergency spillway at the west abut- 
ment, through a gully and into Mono Creek. 

The outlet structure was constructed in a cut in foundation soils near the 
maximum section of the dam. It follows the foot of the right hand slope of 
Mono Creek, from the upstream to the downstream toe of the dam. It con- 
sists of a reinforced concrete culvert, 800 feet long, housing a 54" diameter 
steel pressure conduit, controlled at the downstream end by a Howell-Bunger 
valve. The water is admitted into the steel outlet conduit through a reinforced 
concrete inlet structure. At the intersection between core and outlet conduit, 
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the base of the core is located at a minimum depth of 5 feet below the base of 
the outlet conduit. Therefore at that location the ovtlet conduit is entirely 
surrounded with impervious core material. 


Construction of the Dam 


The construction of the dam was started in May 1953. It required the fol- 
lowing quantities of fill material: 


1,236,824 cu. yds. Core material and blanket under dam 

1,027,208 cu. yds. Upstream blanket material and cutoff 

2,538,319 cu. yds. Semi-pervious and pervious materials for the 
upstream and downstream portion of the dam 


486,546 cu. yds. Slope and blanket protection materials and filter 
materials 


5,288,897 TOTAL (Compacted Volume) 


Prior to construction, the damsite was occupied by an old, high-altitude 
forest. The trees had diameters up to 48" and the ground was locally strewn 
with big boulders. Therefore, the stripping of the site was a major operation. 
The boulders were blasted and all stones larger than 6 inches were removed 
and stockpiled, to be used as riprap. The trees were either uprooted or cut 
at an elevation of one foot above ground level. All stumps and roots were re- 
moved to a distance of 20 feet beyond the outer boundaries of the dam. After 
this preliminary operation, the area to be occupied by the embankment was 
harrowed or disked to a depth of 6 inches and moistened by sprinkling to bring 
the water content of the processed layer up to optimum. Prior to placing the 
embankment the top layer was cross-compacted by 8 passes of a 50-ton rub- 
ber tired roller. The fill material was spread in 8 inch layers and compacted 
in the same manner as the aforementioned top layer. 

Construction proper started in June 1953 by making the excavations for the 
toe drain, the outlet conduit, and for the outlet channel. In order to excavate 
for the outlet conduit, located in the normal course of the stream, it was ne- 
cessary to divert Mono Creek at a point 1200 feet upstream from the upstream 
toe of the dam. The creek was diverted into an open cut channel excavated 
along the foot of the left hand slope of the riverbed. Since the bottom of the 
trench for the outlet conduit, at the site of its intersection with the core of the 
dam, was located about 5' below the water table, the trench was drained by 
pumping from well points. 

In September 1953, the outlet conduit was concreted and the intake channel 
was provided with an “impervious” blanket. At that stage, Mono Creek was 
diverted from its artificial bed into the lined inlet channel and through a by- 
pass pipe into the outlet conduit, thereupon the balance of the blanketing opera- 
tions in the Mono Creek Valley upstream from the dam was completed. Hence- 
forth, Mono Creek flowed through the outlet conduit into the outlet channel 
(see Figure 5), and entered the old stream at a point about 800 feet down- 
stream from the downstream toe. None of its waters entered the toe trench. 

In August 1953, the drilling of the relief wells was commenced. Figure 7 
shows a typical relicf well installation. In drilling the relief wells, it was 
noted that silts from higher strata would cave into the hole. To prevent this, 
a seal of bentonite, 6 to 8 feet thicker than the silt layer was placed around 
the well casing above the pervious strata. This confined the flowing water to 
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the casing, and prevented the fine sands and silts from washing down and into 
the pervious strata which is essential to the efficient operation of the well. 
The flow of water towards the pervious filter well produces a natural filter 
surrounding this section. 


Modifications of Design During Construction 


The excavation for the cutoff was started October 1953 at the upper edge 
of the eastern (left hand) slope of the creek channel. As excavation proceeded 
in an easterly direction, the depth of excavation to reach the silt stratum was 
generally about 22' below the broad flat floor of Vermilion Valley. At a dis- 
tance of about 1400 feet from the east abutment, contact with the silt stratum 
was lost. In order to get reliable information concerning this unanticipated 
feature of the subsurface geology, subsoil exploration by borings was resumed 
in the Spring of 1954. It involved the drilling of 28 boreholes with depths up 
to 90 feet. In many holes, undisturbed samples were taken at points spaced 

5 to 10 feet vertically. On each sample, mechanical analysis and permeability 
tests were performed. 

Some of the 28 holes were drilled along the centerline of the proposed cut- 
off, east of the point where contact with the silt stratum was lost, and the 
others north of the centerline. The results of the borings along or close to 
the centerline are shown in Figure 9. The diagram also reveals the erratic 
character of the pattern of stratification of the fluvio-glacial deposits. The 
results of the borings, combined with those of the preceding geological survey, 
led to the conclusion that the cutoff, for a distance of about 300 feet, is located 
above a gap in the silt stratum, A in Figure 5. This gap was eroded by a melt- 
water stream and then backfilled with fluvio-glacial sediments before the 
youngest till sheet was deposited. During the last advance of the ice, the 
stream channel deposit A was buried beneath a till sheet and a set of terminal 
moraines. 

After the ice retreated from the terminal moraine M in Figure 5, a stream 
channel B was eroded along the inner or northern edge of the moraine and 
backfilled with fluvio-glacial sediments. At a distance of about 600 feet north 
of the cutoff, the bottom of the younger stream channel deposit B intersects 
the top surface of the older one, A. Hence at that point, water can flow out of 
the reservoir through the deposit B into and through the deposit A under the 
dam towards the left hand slope of Mono Creek downstream from the dam. 

Because of this geological situation, it was decided to excavate the cutoff 
trench to a maximum depth of 25 feet along the entire line indicated on Figure 
5 in spite of the fact that over a length of 300 feet the base of the cutoff is 
located on fluvio-glacial sediments and not on silt. However, in order to re- 
duce the quantity of water which will escape out of the reservoir into the 
stream channel deposit A, the section of the cutoff located above the gap was 
supplemented by a blanket extending from the cutoff in an upstream direction 
along B, Figure 5, to a point beyond the area where the stream channel de- 
posits A, and B intersect. The resulting layout for cutoff and blanket is shown 
in Figure 8c. 

Another modification of the original project became necessary when it was 
realized that the available quantity of highly pervious sand and gravel was 
smaller than anticipated. The cost of screening the silty embankment materi- 
al would have been prohibitive. Therefore, it was decided to subdivide the 
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downstream section of the dam into a very pervious and a less pervious sec- 
tion. The very pervious material was placed on the base of the dam and 
against the lower part of the downstream face of the core of the dam as shown 
in Figures 6a and 6b Zone I. Since the estimated line of seepage is located 
entirely within the very pervious material, this modification of the design had _ 
no effect on the stability conditions. 
Finally, the necessity of some minor modifications of the original design 
was realized when the observations during the first filling of the reservoir 
began to reveal the pattern of subsurface drainage. These modifications in- 


volved chiefly the installation of drains in the natural ground downstream 
from the dam. 


Means for Observation 


The performance of the dam has been observed by means of piezometers, 
discharge gages, and settlement reference points. The pattern of subsurface 
drainage and the performance of the relief wells is disclosed by the piezo- 
meter readings. The relationship between the water level in the reservoir 
and the amount of leakage out of the reservoir was determined by correlating 
the measured discharge at the outlets of the drains with the corresponding 
water level in the reservoir. The settlement of the different parts of the dam 
was measured by repeating from time to time the survey of the settlement 
reference points. 

The design of the piezometers had to be adapted to the permeability of the 
soil surrounding the perforated tip and the location of the upper end of the 
piezometer with reference to that of the tip. Therefore, three different types 
of piezometers, A to C, have been used. 


Type A: Observation wells, for measuring the pore water pressures in 
pervious material. They consist of 1 inch diameter tubes terminating in 
Stang well points, installed in 8-in. drill holes. Above the well points, the 
space between tubes and natural ground was filled with compacted clay. 


Type B: Casagrande piezometers, for measuring pore pressures at multi- 


ple elevations in one well within the body of the dam and in the foundation 
soils. 


Type C: Bourdon gages, attached to the outlet ends of tubes extending 
through the side walls of the outlet conduit into the surrounding fill material. 


The initial piezometer installation consisted of 24 piezometers of Type A 
located midway between the relief wells; 7 triple piezometers of Type B with 
tips located below the base of the dam; 5 piezometers of Type B in the subsoil 
downstream from the dam between the outlet channel and the abandoned chan- 
nel of Mono Creek; 4 piezometers of Type B with tips in the core; 8 Type A 
piezometers with tips in the downstream portion of the dam; and 5 Type C 
piezometers in the side walls of the outlet conduit. The total number of 
piezometers is 53. Their locations are shown on Figure 10. 

Each of the 7 triple piezometers of Type B consists of three 1/2" O.D. 
plastic tubes installed in one 4-in. drill hole. The tips are located at about 
elevation 7420, 7450 and 7480. Between tips, the space between tubes and 
natural ground was sealed with clay. The uppermost tips at elevation 7480, 
are located about 15 feet below the level of the original bottom of the valley of 
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Mono Creek at its intersection with the centerline of the dam. During the 
first filling of the reservoir supplementary piezometers were installed at 
those few points where additional information concerning the pore pressure 
conditions was called for. 

When the filling of the reservoir was started, the loss of water due to 
seepage was measured at two points. One of them was located in the short 
section of abandoned streambed of Mono Creek, about 310 feet downstream 
from the lowest point of the toe drain at Station 28. The second one is located 
at about Station 8 at the bottom of a gentle sag in the toe drain between the 
west abutment of the dam and about Station 10. The station 8 toe drain outlet 
consists of a 90° V-notch weir equipped with a staff gage. The outlet at Sta- 
tion 28 consists of an 18" rectangular weir plate and is equipped with an 
automatic recorder. 

During the first filling of the reservoir, water seeped out of the ground at 
three points downstream from the dam. At these points the seepage was col- 
lected in drains and gauging devices were installed at the outlets of the 
drains. 

In addition to the piezometers and gauging stations, settlement reference 
points were established. These are spaced at intervals of 300 feet and con- 
sist of brass markers inserted into the tops of cylindrical concrete blocks 
with a diameter of 12 inches. The base of each block is located 4 feet below 
the fill surface. Ten such points are on the crest of the dam and five on the 
downstream berm. Provisions have also been made for measuring the settle- 
ment of the invert of the outlet conduit. 


Pattern of Seepage 


As was expected on the basis of the geological history of the glacial sub- 
strata of the damsite, the piezometer readings disclosed an intricate, three- 
dimensional pattern of seepage. In some portions, the flow lines rise and in 
others they descend in the direction of the flow. Therefore, at different ele- 
vations the curves of equal piezometric pressure have a different shape. 

This is shown in Figure 10 representing the curves of equal piezometric heads 
at elevations 7420 and 7480 respectively for full reservoir. If the relief 

wells had not been drilled, the piezometric levels at elevations 7420, down- 
stream from the highest portion of the dam, could have been very much higher 
than those shown in Figure 10, involving the possibility of the formation of 
boils on the valley floor, or even of a heave of the floor. Therefore, the in- 
stallation of the relief wells was an essential precaution. 

As the water level in the reservoir rises the piezometric levels in the ob- 
servation wells rise approximately in simple proportion to the corresponding 
rise of the free water level. This relationship is illustrated by Figure 11 in 
which the abscissa represents the elevation, H, of the water level in the 
reservoir and the ordinates of the dash curves the corresponding piezometric 
levels, h, in representative observation wells. Under the following heading, 
it will be shown that the deviations from the straight-line relationship be- 
tween h and H are due to changes in the hydrological conditions beyond the 
area occupied by the reservoir, such as heavy rainfalls or the melting of a 
snow cover, and not due to an erratic scattering from the fundamental rela- 
tionship. 

The slope of the h - H lines in Figure 11, is determined by the value of the 
ratio of the loss of head along the flow line through the piezometric tip 
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between the upper end of the line and the tip, divided by the total loss of head 
along the flow line. With increasing value of this ratio, the slope of the h - H 
line decreases. The slope of the h - H line for the downstream end of any 
flow line is zero. 

Downstream from the core of the highest portion of the dam, the water 
table in the downstream section of the dam is located at an elevation of about 
45 feet above the original ground surface (bottom of the valley of Mono Creek 
prior to the construction of the dam). The slope of the water table, disclosed 
by the piezometer readings indicates that the quantity of water percolating 
through the core is very small compared to the quantity entering the down- 
stream section of this portion of the dam through the bottom and the side- 
slopes of Mono Creek. 

The quantity of water entering the “interior drain” (see Figure 5 and 6b is 
negligible; this fact indicates that the drain serves no useful function. This 
however, could not be predicted. Therefore, the cost of constructing the in- 
terior drain, like that of installing the relief wells adjacent to the two abut- 
ments of the dam, is part of the price which had to be paid for the fact that the 
dam had to be designed and built before the pattern of subsurface drainage 
was known. However, this price was very moderate. 

On the basis of the results of the piezometric observations, the stability 
computations of the downstream portions of the dam were repeated. The 
computations led to the conclusion that the factor of safety Gg of the down- 
stream slope of the dam is consistently greater than the design value Gg = 1.5, 
except for a short section in the proximity of Station 28, where the value of Gs 
might have been slightly smaller than 1.5. Within this section, a small 
amount of fill was added to the downstream section. 


Seepage Losses 


The major portion of the water seeping through the foundation enters the 
toe drain. Originally the toe drain was equipped with two outlets, one at the 
lowest point of the drain, at about Station 28 and another one at about Station 8 
at the foot of a short countergrade in the alignment of the drain. During the 
first filling of the reservoir, the toe drain became obstructed at about Station 
32, probably by an accumulation of iron oxide (FegO3) which is apparently 
leached from the natural soils in the vicinity. This rust-colored, slimy sedi- 
ment collects on the weir plates and in the stilling basins of the outlet weirs 
and therefore is believed to be the cause of the obstruction in the toe drain. 
Most of the water entering the toe drain east of the obstruction comes out of 
the natural ground adjacent to the point of obstruction, Area III, Figure 5. 
There the seepage water is collected in a drain leading to a gauging device. 
Hence, at the present time, the discharge from the toe drain is measured at 
three points. The outlet at Station 8 discharges the water entering the drain 
between Stations 0 and 10 (western section W) the outlet at Station 28, the 
water entering between Station 10 and 32 (middle Section M) and the uninten- 
tional outlet at Station 32 discharges the water entering between Station 32 and 
the east abutment at Station 43 (eastern Section E). 

At full reservoir the discharge measurements at the three gaging points 
furnish the following information concerning the leakage Q into the three sec- 
tions of the toe drain: 
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The Western section, W, 1000 feet long has a hydraulic head of 0 to 60 feet, 
Qw = 0.5 cfs. The unusual features of this section will be discussed below. 

The Middle section, M, 2200 feet long has a hydraulic head of 60 to 142 
feet, Qm = 3.8 cfs. Most of this water comes out of the floor and the side 
Slope of the old stream channel. The contribution of the relief wells to this 
discharge amounts to 0.3+ cfs. 

The Eastern section, E, 1100 feet long, has a hydraulic head of 112 to 0 
feet, Qe = 0.7 cfs. In this section the downstream portion of the dam rests on 
a gentle morainal ridge located between the core and the downstream toe. 
Most of the seepage water flows between the ridge and enters the toe drain a 
short distance above the point of obstruction. 

As the water level in the reservoir rises the rate of flow along the flow 
lines increases in simple proportion to the rise, because the hydraulic head 
increases in simple proportion to the rise of the water level. Furthermore, 
the area flooded by the reservoir also increases; the increase of this area in- 
volves an increase of the area through which water percolates out of the 
reservoir through the subsoil of the dam towards the drainage facilities on the 
downstream side of the dam. Therefore, the loss of water Q out of the reser- 
voir should increase approximately in direct proportion to the square of the 
rise of the water level H - Hy measured from the level of the bottom. The 
Slope of the curve representing the relationship between Q and H - Hg should 
be a curve, the slope of which steadily increases with increasing values of 
H - Ho. In Figure 11, the ordinates of the heavy plain curve represent the 
measured total discharge Q of the toe drain at different elevations H of the 
water level in the reservoir. According to the aforementioned simple rela- 
tionship, the H-Q curve should be a smooth, roughly parabolic curve like the 
dash-dot ideal curve in Figure 11. The height of the shaded areas above the 
smooth curve represents the excess of the measured discharge over the cor- 
responding discharge indicated by the ordinates of the smooth curve. 

An investigation of the sources of the excess discharge indicated by the 
shaded areas led to the conclusion that most of the excess can be accounted 
for by a temporary increase of the quantity of water entering the ground adja- 
cent to the toe drain and beyond the area occupied by the reservoir. The in- 
crease is caused by meteorological events such as heavy rainfall or snow 
melt. As a result of such events, water also enters the downstream portion 
of the dam through the downstream slope, whence it gets into the toe drain. 

In addition to the leakage into the toe drain, seepage also comes out of the 
ground at four places downstream from the dam. Seepage Area I is located in 
the corner between the eastern portion of the dam and the concrete -lined 
spillway channel; Area II is on the bottom of a shallow gully descending to- 
wards Mono Creek from the east, about 500 feet downstream from the toe 
drain; Area III is at the downstream toe near the obstruction in the toe drain; 
and Area IV is in the corner between the emergency spillway and the wester- 
most portion of the dam. The location of these places is indicated in Figure 5 
by the shaded Areas I, II, II, and IV. 

Seepage from Area I, west of the spillway channel did not start until the 
water level in the reservoir approached its highest position. The wet area 
was dried up by the installation of perforated drain pipes, forming a fish bone 
pattern. At the lower end of the drain, a V-notch weir was installed. At full 
reservoir, the maximum discharge is 0.2 cfs. Heavy rainstorms and thaws 
occurring at times of low reservoir stage increase the discharge as much as 
50%. The foregoing facts indicate that the seepage through Area I depends not 


ASCE 


VERMILION DAM 1728-27 


only on the elevation of the water level in the reservoir but also on the eleva- 
tion of the water table beneath the sloping natural ground downstream from 
the eastern portion of the dam. 

Seepage Area II centers about the point where a small spring flowed out of 
the bottom of the shallow gully referred to above. Within this area seepage 
from the reservoir was anticipated. Therefore a stock pile of filter materi- 
als was established in the proximity before the construction of the dam was 
completed. As soon as the water level in the reservoir went up, the discharge 
of the spring increased. Considering the trend of the lines of equal piezome- 
tric level shown in Figure 10, the increase of the discharge can only be ac- 
counted for by the existence of a layer of exceptionally pervious sediments 
leading from the subsoil of the reservoir into the subsoil of the natural ground 
east of the spring. The buried stream channel deposit A in Figure 5 repre- 
sents such a belt. Hence, the increase of the spring discharge accompanying 
the rise of reservoir level confirmed the conclusion derived from the boring 
records, shown in Figure 9, that the belt-shaped deposit A in Figure 5 exists. 

As soon as the discharge of the spring started to increase during the first 
filling of the reservoir, a deep drain was installed on the bottom of the gully 
and the outlet end of the drainage conduit was equipped with a V-notch weir. 
At full reservoir, the maximum discharge of the drain is 0.5 to 0.7 cfs. Dur- 
ing the winter months when the reservoir is at low stage, any unusually wet 
period or snow melt will cause the discharge temporarily to increase by as 
much as 10%. Furthermore, whenever water flows over the service spillway, 
the discharge also increases, indicating seepage from the unlined portion of 
the spillway channel towards Area II. Hence, at times a considerable portion 
of the discharge from Area II can come from sources other than the reservoir. 

Seepage from Area III was discussed earlier in this section. The area re- 
ceives most or all of the water which enters the toe drain between Sta. 32+00 
and 43+00. The leakage from this area started as soon as the quantity of wa- 
ter flowing through the toe drain exceeded the amount which could be freely 
passed through the obstruction located at Sta. 32+00. The seepage from 
Area III, representing an overflow from the toe drain, is in simple proportion 
to the flow through the toe drain easterly from Sta. 32+00 and hence propor- 
tional to changes in reservoir level. 

Seepage from Area IV was not appreciable until the water level in the 
reservoir rose above elevation 7634. However, a further rise of the water 
level was accompanied by a rapid increase of the discharge. In order to 
understand this peculiar phenomenon, the characteristics of the discharge 
from the outlet of the toe drain at Sta. 8+00 should be considered. This outlet 
is located at el. 7572. Yet the discharge from this outlet remained constant 
and equal to about 0.04 cfs. until the water level in the reservoir arrived at 
about el. 7610, which is the elevation of the base of the kame deposit underly- 
ing the west end of the dam (see Figures 1 and 2). Further rise of the water 
level was associated with rapid increase of the discharge at Sta. 8+00. This 
is shown by Curve Cg in Figure 12, representing the relation between the 
water level in the reservoir and the discharge at Sta. 8+00. As soon as the 
reservoir water level rose to about el. 7634, an appreciable amount of water 
also started to emerge from Area IV, and the discharge from this area in- 
creased rapidly with increasing elevation of the reservoir level. The rela- 
tionship between the water level in the reservoir and the discharge from 
Area IV is shown by Curve Cg in Figure 12. The discharge from Area IV is 
represented by the shaded area between curves C2 and C9+C@3 in Figure 12. 
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Both the seepage from Area IV and the abnormal characteristics of the dis- 
charge of the toe drain at Sta. 8+00 can adequately be explained by assuming 
that the permeability of the kame deposit is very much higher than that of the 
fluvio-glacial deposits underlying the eastern part of the dam, including the 
section located between Sta. 32 and 43. Once the water level in the reservoir 
rises above el. 7634 the toe drain no longer suffices to unwater the kame de- 
posit and, as a consequence, part of the water comes out of the adjacent 
ground downstream from the toe drain. While the water level in the reser- 
voir rose from 7634 to its maximum level, 7642, during the first filling of the 
reservoir the discharge increased from 0 to 0.33 cfs. 

The seepage from Area IV was brought under control by the installation of 
an inverted filter covering the area, supplemented by perforated drain pipes. 
The discharge is measured at the downstream end of the collector drain. 

The following is a summary of the results of the discharge measurements 
made during the second filling of the reservoir prior to which time discharge 
from all seepage areas had been intercepted by collector drains and perman- 
ent gaging stations installed. 


Let Q = discharge at full reservoir. 


Toe drain, section W, Station 0-10 Q = 0.55 
Toe drain, section M, Station 10-32 Q= 3.17 
Toe drain, section E, Station 32-43 Q = 0.70 
Seepage Area I (east abutment) Q = 0.17 
Seepage Area II (gully drain) Q = 0.53 
Seepage Area III (in section E above) - 

Seepage Area IV (west abutment) Q = 0.39 


TOTAL Q=5.51 cfs. 


Before the water level reached its highest elevation for the first time, the 
fill material exposed on the downstream slope of the dam became moist along 
several horizontal lines between the highest portion of the dam and the east 
abutment, and at a few points in these lines minor sloughing over areas of 2 
to 3 square feet occurred, but the quantity of seepage responsible for the 
sloughing was imperceptible. The bands were located in the random material 
and generally above the upper boundary of the highly pervious Zone 1 portion 
of the downstream section of the dam (see Figure 6A and 6B). Subsequent in- 
vestigations showed that each one of the moist bands was located above an ex- 
ceptionally silty layer of fill material. On account of the variable silt content 
of the borrow pit material, the formation of such layers was inevitable. The 
situation was easily brought under control by drilling vertical holes through 
the layers down to the pervious Zone 1 fill material and filling the holes with 
clean gravel; thereupon sloughing stopped. 


Settlement 


Settlements are very slight because the site was glacially preloaded, and 
the fill well compacted; and also because the well graded materials of the sub- 
soil and the dam contain little or no clay. 

The maximum settlement of the invert of the outlet conduit is 0.27 foot; and 


since the time when the dam was completed, the crest has settled a maximum 
of 0.12 foot. 
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CONC LUSIONS 


1. If the pattern of stratification of the subsoil of an earth dam is as erra- 
tic as that of the subsoil of the Vermilion Dam information concerning the 
essential features of the subsoil conditions can only be obtained by a thorough 
geological survey supplemented by borings at points to be chosen on the basis 
of geological considerations and by soil tests on representative samples. The 
few borings which were made at the site of the Vermilion Dam fully served 
their purpose. Exploratory work in excess of what has been performed at 
that site would not have furnished any vital additional information to the 
knowledge of the site prior to construction. 


2. The most complete information which can be obtained at a reasonable 
cost concerning significant features of a damsite on a subsoil with an erratic 
pattern of stratification still leaves a wide margin for interpretation. There- 
fore, the construction drawings accompanying the invitation for tenders are 
inevitably based on more or less arbitary assumptions concerning the charac- 
ter of the subsoil. While these assumptions are being made, the engineer in 
charge of design should consider, on the basis of the results of the geological 
survey, all the uncertainties involved, and he should make all the provisions 
which are necessary for detecting the errors in these assumptions before it 
is too late to modify the design during the construction in accordance with the 
findings. The procedure is illustrated by the contents of this paper. 


3. Because of following the procedure described in this paper, the 
Vermilion Dam is as safe and almost as economical as it would have turned 
out to be, if al! the significant features of the subsoil, including the pattern of 
seepage, had been known in advance of construction. The cost of the installa- 
tion of those parts of the dam, such as the interior drain, which were after- 
wards found to be superfluous, were nominal. After the reservoir was filled 
for the first time, there was no gap in our knowledge of the performance of 
the dam, and every phase of its performance is clearly understood. This 
gives the owners and the State Engineer the positive assurance that the dam 
is satisfactory in every respect. 
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PROCEDURE FOR RAPID CONSOLIDATION TEST 
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ABSTRACT 


A rapid method for the consolidation test is suggested. A test of 10 load 
increments can be completed in a few hours instead of 10 days required by 
the conventional method. Precision of results is comparable to that from a 
conventional test and can be checked quantitatively. 


INTRODUCTION 


The standard procedure generally adopted for the consolidation test re- 
quires that each increment of load applied to the test sample is to act for a 
time (usually 24 hours) sufficient to allow the determination of the point of 
100% completion of primary consolidation. Such a long period is inconvenient. 

An alternative procedure for the consolidation test will be discussed in 
this article. The apparatus used for this procedure is the same as that for 
the conventional test. The time required to complete a test series of 10 load 
increments will not exceed one day. 

The theoretical basis of the proposed procedure will be discussed in Part I 


and a suitable testing method is suggested in Part Il. Some illustrations are 
included in Part II. 


Notation 


The nomenclature used in this article is in close agreement with that 
recommended in “Soil Mechanics Nomenclature” published by the ASCE in 
1941. The only exceptions are as follows: 


a,b - dial reading in A-type, B-type test respectively; 
{ - conformity factor; 


Note: Discussion open until January 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1729 is 
part of the copyrighted Journal of the Soil Mechanics Division, Proceedings of the 
American Society of Civil Engineers, Vol. 84, No. SM 3, August, 1958. 
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h_ - difference of 2 intercepts, see Fig. 2; 


- the slope of the steepest tangent on the time compression curve. 


Part I 


Theoretical Basis of the Steepest Tangent Fitting Method 


The equation governing primary consolidation based on Terzaghi’s one- 
dimensional theory is 


(1) 


The solution of equation (1) can be found by means of Fourier’s transform, 
and the average consolidation ratio can then be given as follows 


= = M 


where M = 


5 (m + 1), m is an integer, and € is the base of natural logarithms. 


-2 
Assuming 4= 4 + & a +l Sin fe where uj, Ug and ug are constants (See 


Fig. 1.), we have 


U= 1-5 
meo 


where ¢(M) is a function of M only. 
The residual of consolidation R is 


R=i/-U = 5 


The general expression for the slope of a tangent on the R-logT curve is 


d(/agT) 


lag dT 


yar em (5) 


The condition for the point of inflexion, through which the steepest tangent 
passes, is that 


6) 
d (lagT}* 
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CONSOLIDATION TEST 


C1- MT) = 0 
mao 


From eq. (7), we obtain T = 0.404 which corresponds to R = 0.298 or U = 
0.702. Substitution of this T-value into eq. (5) with P(m)=%2 ie., u7=0 
gives 


S= dR = 0688 
| 720-404 


The values of T and R are practically independent of ¢(M). But s can 
u u 
theoretically vary from 0.650 (with ug = - = ) to 0.715 (with ug = + | a 


Matlock and Dawson gave s = 0.68, which presumably was obtained on the 
assumption of initial linear distribution of pore water pressure. (See Fig. 1 
and Ref. 1) A detailed discussion of the s-value is beyond the scope of this 
paper. Its actual variation fortunately is within a smaller range than that 
mentioned above. The value of 0.688 may be taken as a standard value, al- 


though the use of a smaller and simpler value gs = : is preferable. 


The dial reading, ay, (or void ratio), corresponding to x% consolidation can 
be estimated by means of the following equation 


ial h x (8) 
aA = 100 


where ao is the corrected initial reading, s is the steepest slope of the time 
curve, and h can be obtained from the experimental results. The determina- 
tion of h is illustrated on Fig. 2. 


Log p vs e Curve for Partially Consolidated Soil 


Suppose two different series of oedometer test are performed on a soil 
starting from identical conditions. 

In series A, every load increment is applied when the consolidation process 
is only 90% completed under the previous load. In series B, the load incre- 


ment is applied at the stage of 150% consolidation, which would be the 24-hour 
reading for a soft clay. (See footnote) 


FOOTNOTE: 


The orthodox division of consolidation into the primary and the secondary 
part is rather arbitrary and only suitable for Terzaghi’s rheological model of 
soils. If a rheological model which can fully represent the real soil is found 
and adopted, there will be no “secondary” consolidation at all. Therefore, it 
is preferable to refer to the “secondary” consolidation in terms of a percent- 
age of the “primary” part, which is defined by a mathematical theory. 

The 150% consolidation therefore means a consolidation with a compression 


ratio r = 100/150 = 0.67, assuming that the corrected zero point coincides with 
the initial dial reading. 
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The curves shown in Fig. 3 have been constructed as follows: 


Ao, By are curves representing the initial dial gauge (or void ratio) read- 


ings, ay bos respectively, where the superscript denotes the i-th increment 


of loa-d. Agg, Bj509 are curves representing the final readings, 299) bi go 


respectively. 


These two series of tests were started under the same initial conditions. 


That is, a’ = by = 0 at pressure Py If the samples are identical and consoli- 


date at the same time rates, then the following relation also exists 


(9) 


/ 
Faq = 


The second increment of load is applied in both cases. In series A, 


A . In Series B, be = . Obviously > ,and 62>a? 


(a and b increase with decrease in sample thickness). And in general, 


(10) 
65 24; 


Relation (10) can be replaced by a more useful expression 


(11) 
16 = - 6 


where the constant 1.6 is the ratio of 150% and 90%. The value of - can 


therefore be estimated from the results of test series A. As a consequence, 
the test series B which would have been performed to obtain a curve B 


150’ 
may be replaced by the test series A. 


Part II 


Suggested Procedure for Consolidation Test 


In the following paragraphs, only those steps which depart from the con- 
ventional procedure are discussed. The major difference appear in the stage 


from application of load to the specimen until the plotting of load-compression 
curves. 


The following recommendations are made: 


1. Preparation for plotting dial readings directly on semi-log graph paper 
should be completed before the first load increment is applied. The size of 


the graph paper should be large enough for a precise construction of the steep- 
est tangent. 


2. Dial readings should be taken at the following time intervals: 


every 5 sec. for the first minute, 
every 30 sec. for the next 9 minutes, and then 
every minute till t = to or a= ago whichever occurs later. 


= 
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3. The corrected zero reading a, is determined by the parabolic method (See 
Ref. 2) while the test is proceeding. When the steepest tangent to the time- 
compression curve has been found, the difference in dial reading for one log- 
cycle, h, is determined. The dial reading corresponding to x% completion of 
consolidation is given by eq. (8). (In the following paragraphs whenever dial 
readings are mentioned, they refer to the change in dial reading from ap.) 


4. The last reading should be taken at time t = too unless agg has not yet been 


reached. The value for too can be estimated during the test by the following 


equation: 
(12) 


where teo can be located on the semi-log graph after finding a_,. from eq. (8). 


50 
There is, however, another requirement to be satisfied when taking the 
last reading. This is that a must equal to agg where Ag can be found from 


eq. (8). Theoretically, the point with an abscissa of log t,,, should coincide 


99 
with the point whose ordinate is Ago: If they do not coincide, another incre- 


ment of load should not be applied until both these points have been recorded. 


5. The conformity factor f should then be found. f,,, is defined as the ratio of 


90 
the actual dial reading to the predicted dial reading at 90% consolidation. 

The actual dial reading is found on the time-compression curve by setting 
t= too which is determined from eq. (12). The predicted dial reading is given 


by eq. (8). 


6. The above procedure should be repeated for the following increments of 
load. For the last increment, the load may be left on to obtain the 24-hour 
“final” reading as in the conventional method. 


7. Two laboratory curves Ap and Ago can then be drawn. Other curves can be 


interpolated or exterpolated when required. 


Comments on the Suggested Procedure 


1. Advance preparation for plotting of the dial readings and for the construc- 
tion of the steepest tangent as well as the subsequent calculation is necessary 


because the person performing the test will be fully occupied while the test is 
in progress. 


2. Sufficient readings should be taken to provide a curve enabling a well de- 
fined tangent to be drawn. 


3. The steepest tangent can be constructed when the path of the points on the 
semi-log graph paper indicates a slight decrease in slope. 


4. The point of 90% consolidation is chosen as the last point only for the sake 
of convenience. The point of inflexion appears at about 70% consolidation, and 
the test can be stopped at any time after this. However, all the points around 
the 70% point appear to lie on a straight line, which extends over a range 


J 
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equivalent to nearly 30% of the “primary” compression. The only means of 
ensuring that the steepest tangent point has been passed is to wait for the 
curve to turn the other way. The 80% point may still be in the lower reach of 
the straight portion on a small scale graph and the use of the 90% point is 
very much on the safe side. 

The other advantage of using the 90% point is that the square root of time 
fitting method (See Ref. 3) can then be applied to the test results and a valu- 
able comparison can be obtained. This is, however, not a necessity for a 
routine laboratory test. 


5. The conformity factor is a measure of the discrepancy between the actual 
behaviour of the sample and its expected theoretical behaviour. Theoretically, 
this factor is constant for identical samples under identical conditions, and its 
variation in practice shows the error of observation. The mean square error 
can be found, thus giving a quantitative indication of the precision and uni- 
formity of the experiment. 


6. The 24-hour final reading provides further information on the conformity 
between behaviour of the test sample and the prediction made from the theory. 


Ne Ao and Ago are two basic curves, from which other curves can be con- 


structed. 


Part Il 


Illustration of Steepest Tangent Fitting Method 


A typical oedometer curve of dial reading vs time on a semi-log plot is 
shown in Fig. 4. The sample used was Wyoming bentonite with calcium as its 
exchangeable ion. The liquid limit was 240% and the coefficient of permeabili- 
ty was of the order of 1 micron per second, nearly 1000 times less than that 
of an ordinary clay. 

It can be seen on Fig. 4 that a period of not less than five days was neces- 
sary for the determination of the 100% consolidation point by Casagrande’s 
method (See Ref. 2) The point recorded at t = 110 minutes clearly indicates 
the existence of the steepest slope and the “rapid” test was therefore com- 
pleted in less than two hours’ time. 


A. Steepest Tangent Fitting Method: 


i) Corrected zero reading 4d, = 6570 
ii) Intercepts of the steepest tangent for one log cycle, 


h = 6438 (t = 10) - 6175 (t = 100) = 263 
h = 6450 (t=9) -6187(t=90) = 263 Check. 


iii) Predicted dial reading for 90% consolidation 
dog = 6570 - 1.5 x 263 x 0.90 = 6214 
iv) Actual dial reading for 90% consolidation 
a) de = 6570 - 1.5 x 263 x 0.50 = 6372 
b) from Fig. 4, te corresponding to d = 6372 was 17.8 min. 
c) too = 4.3 te = 76.4 min. 
d) the dial reading at t = 76.5 min. was 6213. 


| 
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v) Conformity factor foo = (6570 - 6214)/(6570 - 6213) = 1.003 
B. Casagrande’s Method: 

i) d, = 6570 

ii) dygo = 6161 
C. Comparison of values of dj Qo: 

i) Method A, qi 00 = 6570 - 1.5 x 263 = 6175 

ii) Method B, 00 = 6161 


iii) Discrepancy: (6570 - 6153)/(6570 - 6161) = 1.02 


Illustration of the Use of A-Curves 


Results of an oedometer test performed in accordance with the suggested 
procedure (except that the test was stopped at 80% consolidation instead of 
90%) are shown in Table 1. The corresponding curves Ao and Ago are shown 
in Fig. 5. 

It should be noted that in this example, the initial value of the void ratio 
was equal to the final void ratio under a previous load. This occurred be- 
cause the sample was prepared under a high vacuum, such that there was 
practically no air entrapped in the sample. If there had been air in the sam- 
ple, there would have been instantaneous consolidation from compression of 
the air, which would have had to be included in the calculation. 

The calculation used in development of Table 2 was as follows: 


A. Rearrangement of Results. (Column 1 to Col. 5). 


i) Col. 1 to 3 were taken from Table 1. 
ii) Col. 4 = Col. 2 - Col. 3, ie., (4e)) = @'- en 


Arithmetic check 6, - ef = 5 (Le), + 


1.183 - 0.830 = 0.358 Check. 
iii) Col. 5 = (Col. 4) i.e., = 
358 x 15/8=671 Check. 


B. Preparation of B-Test Data. (Col. 6 & 7). 


i) The same sample was assumed to undergo a B-test. The initial con- 


dition was pi = e! = 1.188 


ii) Col. 7= Col. 5 - Col. 6 ice., - = 


or, + )+ = bi + = 


SO 


if #0 
Arithmetic check: - bf, = 
1.188 - 0.517 = 0.671 


| 
Check. 
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C. Construction of B-Curves 


in Fig. 5. 


The construction of B-curves using the data in Col. 1, 6, and 7 was shown 
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In a case where the original oedometer data are not available, a graphical 
method comparable to the above calculation can be devised. However, it 
would be better to read off from the curves and derive data from them than to 


rely entirely upon the graphical method. 
Table 1 


A-Test Data. p vse. 


Initial Void Ratio 


1.188 
1.184 


Applied Pressure 


wro 
3 


CONCLUDING REMARKS 


which unnecessarily long when the conventional method is used. 


Final Void Ratio 


1. The “rapid” procedure can shorten the duration of a consolidation test, 


1.184 
1.179 
1.173 


1.161 
1.107 
1.037 


0.989 
0.899 
0.820 


2. The reliability of the rapid procedure is not inferior to that of the conven- 
tional procedure. The rapid method can indicate the precision of a series of 
tests quantitatively, whereas the reliability of the results obtained by the con- 


ventional method can only be found qualitatively. 


| 
Table 2 | 
Column No. 1 2 3 4 5 6 yf 
Appl. 
0.1 1.188 1.184 4 8 1.188 1.180 
0.2 1.1799 5 9 1.180 
0.4 1.179 1.173 6 Aa 1.272 1.160 
0.8 1.173 2.561 12 23 1.160 1.137 
1.6 1.161 1.107 54 101 LAS? 1.036 
3.2 1.107 1.037 70 131 1.036 0.906 
6.4 1.037 0.969 68 128 0.905 0.777 
12.8 0.969 0.899 70 131 0.777 0.646 
25.6 06.899 0.820 69 129 0.646 0.517 . 
Sum = 358 671 
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3. The steepest slope fitting method, although shown in this article to be 
based upon the classical rheological model of clay suggested by Terzaghi, can 
be applied to other complex models such as Biot’s or Tan’s generalised 
model, The only differences are in the s-value and in the coordinates of the 
tangent point. (See Ref. 4 and 5 for details of these models.) 


4. After some slight modification, the square root of time fitting method, 
which also makes use of the steepest tangent, can be incorporated in the 
“rapid” method to provide a further check. 
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EFFECTS OF GROUND ON DESTRUCTIVENESS OF LARGE EARTHQUAKES 


C. Martin Duke,! M. ASCE 
(Proc. Paper 1730) 


SYNOPSIS 


The present state of empirical knowledge of the relation between ground 
conditions and earthquake damage to structures is presented in terms of 
observations in destructive earthquakes. Thirty-six earthquakes are cited for 
which a relation has been reported. Patterns of behavior exist for several 
kinds of building. Firm ground is much preferred to soft ground for most 
kinds of structure. Strong rigid buildings may be an exception. 


INTRODUCTION 


One of the major considerations facing the engineer who designs structures 
for seismically active regions is the question of how the ground under a 
structure influences the probability of damage in an earthquake. This question 
has interested the observers of earthquakes for at least three centuries. The 
earthquake-ground-structure system may be visualized as a circuit witha 
feedback loop from structure to ground. 

Scientific investigations have taken three forms: the recording of relations 
between ground conditions and destruction in large earthquakes, field and labo- 
ratory studies using small natural and simulated earthquakes, and theoretical 
work, This paper summarizes the present state of knowledge gained from 
large earthquakes. The several figures presented are idealizations of the re- 
ports of the respective authors, to whose papers the reader is referred for 
details. 

Thirty-six earthquakes for which a relation has been reported between 
damage and ground conditions are listed in Table I. 


Note; Discussion open until January 1, 1959. To extend the closing date one month, 

a written requést must be filed with the Executive Secretary, ASCE. Paper 1730 
is part of the copyrighted Journal of the Soil Mechanics Division, Proceedings of 
the American Society of Civil Engineers, Vol. 84, No. SM 3, August, 1958. 


1. Dept. of Eng., Univ. of California, Los Angeles, Calif. 
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TABLE I 


THIRTY-SIX EARTHQUAKES FOR WHICH A RELATION HAS BEEN REPORTED 


BETWEEN BUILDING DAMAGE AND GROUND CONDITIONS 


Place Bibliographic 
Citation 


Port Royal, Jamaica 2 
Lisbon, Portugal 

Lisbon, Portugal 
Messina, Italy 

Calabria, Italy 
Talacahuano, Chile 

Tokyo, Japan 

New Zealand 

Mendoza, Argentina 

San Francisco, California 
Yokohama, Japan 

Manila, Phillipine Islands 
Charleston, South Carolina 
Mino-Owari (Nobi), Japan 
Assam, India 

Valparaiso, Chile 

San Francisco, California 
Kingston, Jamaica 
Germany 

Freiburg, Germany 

Tokyo, Japan 

Long Beach, California 
Oberschwaben, Germany 
Tottori, Japan 

Tonankai, Japan 

Mikawa Province, Japan 
Nankaido District, Japan 
Fukui, Japan 

Seattle, Washington 
Imaichi, Japan 

Assam, India 

Hokkaido, Japan 

Kern County, California 
Eureka, California 
Quetta, Pakistan 

Oaxaca, Mexico 


1692 
1737 
1755 
1780 
1.783 
1835 
1855 
1855 
1861 
1868 
1880 
1880 
1886 
1891 
1897 
1906 
1906 
1907 
1911 
1913 
1923 
1933 
1935 2 
1943 3 
1944 15 
1945 3 
19446 3, 16 
1948 7, 18 
19h9 9 
1950 
1952 2 
1952 3 
1954 an 
1955 PS 
1957 
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History to 1886 


The following quotation from John Milne, engineer and seismologist, (1) 
summarizes the early qualitative findings. 


“That the nature of the ground on which a building stands is intimate- 
ly related with the severity of the blow it receives is a fact which has 
often been demonstrated . . . 

“In the great Jamaica earthquake of 1692, the portions of Port Royal 
which remained standing were situated on a compact limestone foun- 
dation; whilst those on sand and gravel were destroyed . . . According 
to the observations made at Lisbon, in 1737, by Mr. Sharpe, the destroy- 
ing effects oi this earthquake were confined to the tertiary strata, and 
were most violent on the blue clay, on which the lower part of the city 
is constructed. Not a building on the secondary limestone or on the ba- 
salt was injured. 

“In the great earthquakes of Messina, those portions of the town situ- 
ated on alluvium, near the sea, were destroyed, whilst the high parts of 
the town, on granite, did not suffer so much. Similar observations were 
made in Calabria, when districts consisting of gravel, sand, and clay 
became, by the shaking, almost unrecognizable, whilst the surrounding 
hills of slate and granite were but little altered. At San Francisco, in 
1868, the chief destruction was in the alluvium and made ground. 

“At Talacahuano, in 1835, the only houses which escaped were the 
buildings standing on rocky ground: all those resting on sandy soil were 
destroyed... 

“Humboldt observed that the Cordilleras, composed of gneiss and 
micaslate, and the country immediately at their foot, were more shaken 
than the plains. 

“Some writers have asserted that the wave-like movements (of the 
Calabrian earthquake in 1783) which were propagated through recent 
strata from west to east, became very violent when they reached the 
point of junction with the granite, as if a reaction was produced when the 
undulatory movement of the soft strata was suddenly arrested by the 
more solid rocks. Dolomieu when speaking of this earthquake says the 
usual effect was to disconnect from the sides of the Apennines all those 
masses of sand and clay which either had not sufficient bases for their 
bulk, or which were supported only by lateral adherence... 

“When referring to the question as to whether buildings situated on 
loose materials suffered more or less than those on solid rocks, Mallet, 
in his description of the Neapolitan earthquake of 1857, remarked: ‘We 
have in this earthquake, towns such as Saponara and Viggiano, situated 
upon solid limestone, totally prostrated; and we have others such as 
Montemarro, to a great extent based upon loose lays, totally levelled. 
We have examples of almost complete immunity in places on plains of 
deep clay as that of Viscolione, and in places on solid limestone, like 
Castelluccio, or perched on mountain tops like Petina.’... 

“Professor D. S. Martin, writing on the earthquake of New England 
in 1874, remarks that in Long Island the shock was felt where there was 
gneiss between the drift. . . . generally the shocks were felt more 
strongly and frequently on rocky than on soft ground... 


} 
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“... a civil engineer, writing about the New Zealand earthquake of 
1855, when all the brick buildings in Wellington were overthrown, says 
that it was most violent on the sides of the hills at those places, and 
least so in the centre of the alluvial plains... 

“, . . in the destructive earthquake of 1855, when a large portion of 
Tokio was devastated, it was a fact, remarked by many, that the dis- 
turbance was most severe on the low ground and in the valleys, whilst 
on the hills the shock had been comparatively weak... 

“In Hakodadi, which is a town situated very similarly to Gibraltar, 
partly built on the slope of a high rocky mountain and partly on a level 
plain, from which the mountain rises, the rule is similar to that for 
Tokio, namely, that the low flat ground is shaken more severely than the 
high ground. At Yokohama, sixteen miles southwest from Tokio, the rule 
is reversed, as was very clearly demonstrated by the earthquake of 
February 1880, when almost every house upon the high ground lost its 
chimney, whilst on the low ground there was scarcely any damage done; 
the only places on the low ground which suffered were those near to the 
base of the hills. The evidence as to the relative value of hard ground 
as compared with soft ground, for the foundation of a building, is very 
conflicting. Sometimes the harder ground has proved the better foun- 
dation and sometimes the softer, and the superiority of one over the 
other depends, no doubt, upon a variety of local circumstances .. .” 


Recent History 


John R. Freeman, civil engineer,(27) has presented a good summary of 
worldwide experience, with emphasis on the early 20th century. The Report 
of the State Earthquake Investigation Commission, (5) the reports of the Im- 
perial Earthquake Investigation Commission, (10) and the reports of the 
Japanese committees(17) and the U. S. Army,(18) give extensive information 
on the ground-intensity correlations of the earthquakes of San Francisco 1906, 
Tokyo 1923, and Fukui 1948, respectively. Hisada and Nakagawa(28) present- 
ed a detailed summary and bibliography of Japanese findings since 1945. The 
author(29) has compiled an annotated bibliography of worldwide experience. 
The first detailed quantitative study was made by Wood on the 1906 San 
Francisco earthquake. (5) For most subsequent earthquakes such quantitative 
studies have been made, presenting comparisons of distribution of degree of 
damage to various kinds of buildings with distribution of ground conditions. 
The recording of accelerations in strong earthquakes, initiated in 1932 
in the United States and in 1952 in Japan, has helped clarify the effects of 


ground on intensity. Records have been obtained to date in a number of strong 
western United States earthquakes. 


Japanese Houses 


Extensive investigations of the relation between ground conditions and 
damage to wood houses have been conducted in Japan since 1923. The usual 
criterion of destructiveness was the “damage ratio” or ratio at a given lo- 
cality of damaged houses to total houses. Ground parameters used by various 


investigators have been age of surface deposits, depth of alluvium, softness 
of ground, and seismic wave velocity. 


| | 
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The older Japanese houses, which constitute the majority of structures 
damaged in past earthquakes, are constructed of wood, usually with heavy 
roofs of tile or thatch, and usually with exterior walls of mud plaster over 
bamboo lath. Principal lateral force resistance comes from these walls, the 
timber framing and interior partitions being of such construction that they 
provide only very small lateral resistance. The buildings may be of one or 
two stories. Their average natural period is about 0.3 second, and they have 
been subjected in earthquakes to destructive accelerations of 0.3 to 0.4 of 
gravity. (30) 

Aseismic properties of these houses may be compared with modern United 
States wood-frame construction on the basis of side-by-side performance of 
the two kinds in Fukui in the 1948 earthquake, (18) The latter were designed 
with special attention to earthquake resistance by the Office of the Chief of 
Engineers, Far East Command, U. S. Army, Tokyo. The U. S. style houses 
were damaged in all cases but none of them collapsed or suffered significant 
damage to the structural framing. In the same area of Fukui, the average 
percentages of light-roofed (newer) and heavy-roofed (older) Japanese houses 
collapsed were about 50 and 90, respectively. In the 1935 Formosa earth- 
quake, (31) this kind of Japanese house was much less damaged than were 
adobe houses. 

The great majority of buildings in Japan are located along the seacoast and 
in valleys in the very mountainous topography. The following kinds of foun- 
dations are common: 


1. river deltas, drowned valleys, reclaimed lagoons, muddy alluvium, and 
made land; 

2. littoral sand dunes, beaches, bars, spits, river flood plains, and diluvial 
volcanic formations; 

3. firm tertiary formations, compact conglomerate, and rock. 


The first kind is found most often and the third kind least often. 

It is seen from Table II, where the several recent comparisons of damage 
with ground conditions are summarized, that the three kinds of foundation 
listed above are in the order of increasing desirability. Damage was always 
many times greater on kind 1 than on kind 3, with kind 2 being intermediate 
but usually closer to 3 than to 1 in desirability. 

Certain investigators have made use of numerical ground parameters to 
study the distribution of damage ratio. Fig. 1 shows the influence of depth of 
alluvium found by Kawasumi(32) and Omote and Miyamura.(33,34) The wide 
divergence between the curves for Tokyo and Yokohama may be due to differ- 
ences in quality of the alluvium or to the fact that Yokohama was much closer 
to the epicenter. At Fukui, 1948, (18) damage was roughly proportional to 
thickness of alluvium, with some exceptions. 

Kitazawa(38) obtained a correlation, Fig. 2, between damage ratio and soft- 
ness of ground, He utilized the logs of many borings in which standardized 
blow-count diagrams were obtained, and defined ground softness as the re- 
ciprocal of the area of the blow-count diagram down to a particular depth. 

Omote and Miyamura(34) obtained approximate data on the effect of age of 
made land in Nagoya, Fig. 3. Kanai(20) also found at Imaichi that damage 
diminished with age of reclaimed land, 

Seismic exploration was used by Tanabashi and Ishizaki(13) to find stratifi- 
cation and wave velocities in four cities which experienced destructive earth- 
quakes in 1943-48, After finding no correlations between damage and several 
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TABLE II 


SUMMARY OF EFFECT OF AGE AND TEXTURE OF GROUND ON 


EARTHQUAKE DAMAGE TO JAPANESE HOUSES 


Year Place Ground Damage Ratio, Bibliographic 
percent Citation 
1923 Tokyo (b) Tertiary 0.8-1.)) 32 


Diluvium 0.3-0.5 
Alluvium 223-3 

Yokohama (°) Sand Less than 50 33 
Other alluvium 


1944 Nagoya Ce) Hills 1 3k 


(a) Plain 

Totomi Rock ly 
Sand and gravel 
Clayey alluvium 

Mikawa High ground 5 
Lowlands Greatest 

Tonankai Firm Slight 
Alluvium Severe 
Made land Extra severe 


19146 Kochi Rock or gravel Negligible 
Beach sand Variable 
Swampy ground 

1948 Fukui (4) Hills 


Foothills, beaches 
Alluvial plain 


18 


1952 Hokkaido (») Beach sand 
Peat beds 
Kushiro Shaly hill 
Peat beds 
Made land Greatest 


Damage criteria were: 
(a) Number of collapsed houses. 
(b) Number of collapsed plus number of partially collapsed houses. 
(c) Number of collapsed plus # number of partially collapsed houses, 
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velocity parameters, they obtained a linear correlation, Fig. 4, with the ratio 
of adjusted depth of soft soil to velocity of surface wave. The slope is differ- 
ent for each earthquake, and it would be useful to attempt to reconcile these 
differences in terms of differences in mean intensity and ground quality. 

Suspecting a pattern to the variability between damage rates based on col- 
lapsed and partially collapsed houses, Kanai, (39,40) studied the damage sta- 
tistics of Japanese earthquakes from 1891 to 1947. He found that while the 
great majority of houses on firm ground were not materially damaged, those 
that were damaged were likely to collapse in over half the cases. Conversely, 
damage to houses on soft ground was widespread, but less than half of those 
damaged suffered collapse. He attributes this difference to the fact that differ - 
ential settlement due to the earthquake plays an important role in damage on 
soft ground, while damage on firm ground is primarily due to vibrational 
forces. Such differential settlement has been identified at Hokkaido, 1952, (37) 
and Eureka, 1954, 24) as an important cause of damage on very soft ground. 
This suggests that ground softness, as measured by Kitazawa or by Tanabashi, 
may in some cases be better than depth of alluvium as a criterion of earth- 
quake vulnerability. 


Other Buildings 


Stone and Earth Houses 


Masonry houses of stone, mud, or burned brick, bonded with mud or with 
very weak lime mortar, and houses built wholly of earth, are common in many 


parts of the world. 

The largest amount of information comes from Italy(2,27) from the earth- 
quakes of 1780, 1783, 1860, 1883, 1887, 1894, 1909, and 1930. The typical 
Italian house was two stories in height with walls of stone rubble masonry set 
in weak lime-and-earth mortar, with arches of the same material supporting 
the second floor and with tile roofs on wood framing. In each earthquake list- 
ed, the superiority of rock or other firm ground over soft alluvium was noted. 
Unfortunately, no quantitative studies of this relation have been reported. 

In the earthquakes of 1935 and 1955, the damage in Quetta, Pakistan, was 
greater in the southwestern section on water bearing alluvium than in the 
northeastern section on dry alluvium and gravel.(25) Most houses in 1935 
were of sun dried or burned brick bonded with mud or lime mortar, but those 
erected subsequently were mostly of aseismic design. 

Schlocher and Radbruch(23) made observations in the epicentral area of the 
Kern County, California, earthquake of 1952. Stone masonry and adobe build- 
ings on rock or thin soil were undamaged or only slightly damaged, while 
those on thick soil were damaged severely, damage increasing with thickness 
of alluvium. 

In Chilean earthquakes(41) the behavior of adobe construction has been good 
on rock and gravel ground and very questionable on sandy and clayey soils. 


The best of adobe buildings can not withstand intensities greater than 9 on the 
Sieberg scale. 


Rigid Buildings 


There is evidence that rigid brick and concrete buildings erected since 
about 1850, with strong materials and good workmanship, and in accordance 


| 
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with engineering practice of the time, may be less damaged in earthquakes if 
they are on softer soils. Milne noted(2) that in the New Zealand earthquake, 
1855, brick buildings suffered greatest damage when on hillsides, least when 
on alluvial plains. Also, in the 1880 Yokohama earthquake, almost every 
house on high ground lost its chimney while damage on the low ground was 
very slight. In the 1950 Assam, India, earthquake,(21) brick buildings suffered 
greater damage when on firm than when on soft ground. 

The 1923 Tokyo earthquake has yielded the most quantitative information 
on this subject. In Tokyo there were, on the soft ground, 582 brick and 418 
reinforced concrete buildings, and on the firm ground, 192 brick and 158 rein- 
forced concrete buildings. They were generally of one to three stories. 
Suyehiro(42) reported that the damage to the brick buildings was more severe 
when they were on firm ground. Kanai(43) restudied the data in detail and 
confirmed this finding. In a later study,(44) Kanai showed that the damage 
ratio, based on completely plus partially destroyed brick buildings, was less 
where the thickness of alluvium was greater, as shown in Fig. 5. 

Martel(11) obtained a similar result based on his study of damage to 1261 
one- to four-storied brick buildings in the 1933 Long Beach, California, earth- 
quake, Using the ratio of cost of repairs to original value of building as the 
damage criterion, he determined that damage to buildings on soft water-logged 
soil and on beach sand was somewhat less than to those on more firm ground, 

Damage to the reinforced concrete buildings in Tokyo, 1923, has been ana- 
lyzed by Kanai(43,44) and Kitazawa,(46) Fig. 6. Kanai found that damage di- 
minished with thickness of alluvium, similarly to the case of brick buildings 
but at a lesser rate. Kitazawa plotted his ground softness parameter (see 
above) against damage ratio based on damage to structural elements only. It 
is seen that damage determined on this basis increased with ground softness. 

Japanese experience with small rigid farm structures used for storage 
usually has been consistent with the above. These buildings, called godowns, 
are constructed of closely spaced strong wood pillars covered outside with 
thick mud plaster walls or soft stone and usually have stronger foundations 
than the houses. Small-motion period is about 0.2 second.(47) In the Northern 
Musashi earthquake of 1931, the walls of such structures on paleozoic for- 
mations were cracked to a greater extent than for those on alluvial ground.(42) 
However, in Fukui, 1948, an opposite trend was noted:(48) houses were 
damaged more than godowns when on firm soil, less when on sand dunes, and 
to the same degree when on alluvium. In the 1949 Imaichi earthquake ,(20) 
3174 relatively rigid soft-stone granaries were damaged. While the compara- 
tive damage ratios on various kinds of ground were not reported, it was de- 
termined that the granaries on soft ground suffered most damage in the upper 


1/3 of their height, and those on hard ground were most damaged in the lower 
1/3. 


Buildings in General 


Many other reports have been made of the effect of ground on intensity. 4 
When buildings in general are considered, the superiority of firm ground is 
nearly always proven, This trend appears in most of the reports of the earth- 
quakes of the 18th and 19th centuries, Table I. 

Kitazawa‘46) found for steel frame buildings in Tokyo, 1923, the same 
trend in relationship between damage and ground softness as illustrated in 
Fig. 6 for reinforced concrete buildings. In the Mexico earthquake of 1957,(26) 
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flexible multistory buildings in Mexico City 200 miles from the epicenter 
suffered much greater damage than similar buildings in Acapulco which was 
60 miles from the epicenter. Mexico City is on deep soft alluvium and 
Acapulco is on rock and beach sand. 

Wood’s study(5) of San Francisco, 1906, has become a classic. The San 
Andreas fault was close to the city, and damage had some relation to distance 
from the fault. However, the effect of kind of ground was more pronounced. 
He found the most severe damage on made land, with damage successively 
less severe on alluvium, conglomerate and shale, and rock. Damage was 5 to 
10 times greater on made land than on rock. His findings were supported by 
other investigators. (49,50) 

The same trend has been found for earthquakes in eastern Canada, (51) 
Germany,(12) Chile,(41) Turkey, Greece, and Algeria,(52) and China. (53) 

Neumann(54,55) has provided ground-intensity data on several United States 
earthquakes by drawing envelopes of intensity versus distance from epicenter, 
such as in Fig. 7. He believes that in virtually al) cases the points on the 
lower envelope are for locations on basement rock outcrops, giving a very 
uniform intensity-distance relation which depends only on the character of 
rock (igneous or sedimentary). Points on the upper envelope are believed to 
correspond with locations on soft overburden, with points between the two en- 
velopes representing local geology of intermediate firmness. Neumann has 
found it possible to construct these envelopes from a small number of points 
using curve shapes determined from well-documented earthquakes. 


Strong-Motion Accelerograms 


The strong motion seismograph program of the United States Coast and 
Geodetic Survey has been in operation since 1932, and a large number of 
records has been obtained. In 1956 there were 71 accelerograph stations at 
56 sites in the western United States and in Central and South America, the 
majority being in California. Most are located in the basements of buildings; 
a few are on higher floors, and two are on the ground away from buildings. 

Of the 37 U.S.C.G.S. sites in California, 27 are on unconsolidated sediments, 
4 are on consolidated sediments, and 6 are on rock, However, no adequate 
investigation has been made of the ground conditions at these sites. In ad- 
dition to the U.S.C.G.S. stations, other seismologists operate a few strong 
motion stations in California. 

The Japanese strong motion program was initiated in 1952, but, to date 
there have been only relatively small earthquakes recorded. In 1956 there 
were 14 stations located in seven buildings in Tokyo. The number of instru- 
ments is being increased as rapidly as possible. Although records from the 
present instruments may be influenced by the vibrational properties of the 
buildings housing them, these buildings are situated on various kinds of 
ground, and a large earthquake should yield records with useful data on the 
effect of ground. 

Engineers at California Institute of Technology(56,57) have obtained 88 
spectra based on U.S.C.G.S. records in 11 western United States earthquakes. 
The spectra are of two kinds, velocity and acceleration, and were computed 
for various damping factors including zero damping. 

The damped velocity spectra are consistent with the random disturbance 
hypothesis of earthquake motion(58) in that, above a certain minimum period, 
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velocity is distributed more or less randomly about a mean value of velocity 
which is independent of period. Such a shape suggests the equipartition of 
energy among the earthquake wave components at the surface. Kanai(59) 
shows from amplificatien theory that velocity spectra on very deep alluvium, 
such as exists at most of the stations, should be expected to be flat. 

In several cases, two or more earthquakes were recorded at the same 
station, but examination of the spectra did not disclose any special features 
that could be attributed to the locality. Also, several earthquakes were re- 
corded at two or more stations, but no characteristic resemblance attributa- 
ble to the particular earthquake could be found in the spectra. 

However, in two special cases, features associated with the locality were 
suggested. The station at Helena, Montana, was on rock; the 1935 acceleration 
spectra showed a predominance of short period waves. The station at Seattle, 
Washington, was on watersoaked filled land; the 1949 acceleration spectra 
showed a predominant period of 0.9 second, although the corresponding ac- 
celerations were much smaller than those recorded on firm alluvium at 
Olympia at about the same epicentral distance. 

Damping was introduced so that the spectra might represent the response 
of typical real structures without further computation. However, for analyses 


treating the complete earthquake-ground-structure system, the undamped 
spectra will be needed. 


Lateral Force Coefficients 


Coefficients Based on Earthquake Experience 


Only two cases were found where the variation in equivalent horizontal 
static force with kind of ground has been numerically computed from damage 
studies of particular earthquakes. Kawasumi(61) reconsidered the damage 
data reported by Wood for San Francisco, 1906, to correct for the effect of 
distance from fault. The relative foundation coefficients thus determined are: 


Rock 1 
Alluvium 1.5 
Made land 5 


Takahasi(17) obtained the following relative coefficients for Fukui, 1948: 


Tertiary 0.4 
Diluvium 0.7 
Alluvial ground 1.0 


Marshy land, approximately 1.5 


Coefficients in Building Codes 


In many places the building codes provide for variation in lateral force 
coefficient with kind of ground and kind of foundation. These code provisions 
have been arrived at by committees of engineers who have studied the earth- 
quake experience in their own and other localities. Several examples are 
given below, the coefficients being expressed in fractions of gravity acceler- 
ation. These coefficients, which are to be applied to dead load and a fraction 
of live load, are listed here only to illustrate the present discussion. They 


must not be utilized in design without study of the detailed provisions in the 
referenced publications. 
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Reinforced 
Wood Steel Frame Concrete 


Construction Construction Construction 

Tertiary 12 

Diluvium .16 

Alluvium 5-30 meters deep ‘ .20 
Alluvium over 30 meters 

deep, or made land .20 -20 


-18 
.20 


Seismic zoning maps have been prepared for several Japanese cities, (28) 


Chile(41) Semi -Rigid 
Rigid Construction, Construction, 
Period < 0.4 sec, Period 0.4-0.75 sec. 
Rock .08 .05 
Conglomerated or very compact soil 12 -10 
Sand and filled ground 
With mat or similar foundation -10 .12 
Without mat 


Turkey(52) 


Best Ground Worst Ground 
Higher seismicity zone .02 .04 
Lower seismicity zone .O1 .03 


California Schools (64) 


Safe Bearing Capacity, 
tons /sq. ft. Coefficient 


4 or more .06 
2 to 4 .08 
Less than 2 .10 


CLOSURE 


The type of soil on which a structure is founded, the structure and depth of 
underlying rock, the depth and type of foundation, and the features of surface 
topography including adjacent buildings all have been found to affect the 
damage done to buildings by large earthquakes. While in many of the cases 
cited the soil properties and configuration are quantitatively known, the 
present state of knowledge is such that terms such as “soft ground” must be 
used in generalizations. 

The preponderance of evidence indicates that the majority of structures 
are least damaged when on firm ground. Damage has generally been greater 
on average ground, Soft ground has nearly always been associated with high- 
est damage. 

It must be emphasized that the mass of older damage data were obtained 
on the basis of residential construction, wood in Japan and wood or masonry 
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in the Mediterranean area. There are indications in the more recent data that 
strong rigid structures may be less damaged when on softer ground. 

As noted in the introduction, this paper has considered but one of the three 
avenues of investigation. Such large earthquake studies in the future will be 
more productive if planned with full awareness of the previous work. But the 
experimental and theoretical work going on in this field, principally in Japan 
but also in other countries, is equally important and necessary. There re- 
mains much divergence of opinion among engineers and seismologists as to 
how the kind of ground affects the response of structures, and as to how the 
effects may be quantified. 

The author desires to express his deep appreciation, for assistance and ad- 
vice leading to the preparation of this paper, to Kiyoshi Kanai and Hirosi 
Kawasumi of the Earthquake Research Institute, Tokyo University, to 
Toshihiko Hisada and Kyoji Nakagawa of the Building Research Institute, 


Japanese Government, and to the many other Japanese engineers and seis- 
mologists who aided him. 
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PRESSURE GROUTING FINE FISSURES 


Thomas B. Kennedy! 
(Proc. Paper 1731) 


SYNOPSIS 


Standard field equipment and methods were used in grouting fine fissures 
between specially prepared concrete slabs in a laboratory study with various 
grout mixtures and pressures. The lowest water-cement ratio that could be 
used in grout with neat cement, cement plus fly ash, cement plus fly ash plus 
an intrusion aid, cement plus an intrusion aid, and cement plus calcium ligno- 
sulfonate, to penetrate fissures of 0.01, 0.02, and 0.03 in. width at 25, 50, and 
100 psi was determined. The penetration characteristics and the quality of 
hardened grout films were determined with seven different grout compo- 
sitions, neat cement, cement plus fly ash, cement plus fly ash plus calcium 
lignosulfonate, cement plus calcium lignosulfonate, cement plus finely ground 


water-quenched slag, cement plus pumicite, and cement plus finely ground 
calcined shale. 


INTRODUCTION 


Purpose of Investigation 


It is recognized that the strength and impermeability of a grout film 
follows the water-cement ratio law as do mixtures of other materials con- 
taining water and cement. Because of this fact it is desirable, where im- 
portant work is concerned, to use grouts containing as low water contents per 
unit volume as possible. Since the lower the water content the thicker the 
grout, it is also important to determine the thickest consistencies that will 
penetrate cracks of various widths at given pressures. 


Note: Discussion open until January 1, 1959. To extend the closing date one month, a 
written request must be filed with the Executive Secretary, ASCE. Paper 1731 is 
part of the copyrighted Journal of the Soil Mechanics and Foundations Division, 
Proceedings of the American Society of Civil Engineers, Vol. 84, No. SM 3, August, 
1958. 

1, Chf., Task Committee on Cement Grouting, Concrete Div., U. S. Army 
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In sealing fine seams, grouts having water contents as high as 20 to 1 (20 
cu ft of water per bag of cement) are sometimes used, and grouts of a con- 
sistency as thin as 10 to 1 are not uncommon. The assumption has been that 
thin grout placed under sufficient pressure to force out the excess water used 
to obtain such a consistency will form a hard durable film. Whether it is al- 
ways possible to exert enough pressure to squeeze out this excess water is 
questionable, and for this reason a knowledge of the physical nature of the 
grout films produced by fairly high water-cement ratios would be valuable. 
Such grouts are believed to be extremely pervious and, in time, susceptible to 
leaching that will transform an impervious area into a pervious one. 

The purpose of the program was to obtain information on the degree to 
which the penetration of fine fissures by grout was influenced by surface 
texture of specimen, pumping pressure, water-cement ratio, chemical fluidi- 
fiers, and finely divided mineral additives. It was also desired to determine 
the effect of these factors on the quality of the hardened grout films. 


Description of Investigation 


The study was conducted in three phases during which the ability of grouts, 
of various compositions, to penetrate fine fissures at several pumping 
pressures was observed. Observations were also made of consistency, bleed- 
ing characteristics, and setting time of the grouts. Hardened grout films 
were examined for apparent quality and solubility of selected grout films in 
distilled water was measured. 

The information in the report has been grouped so that all data bearing on 


one condition are presented together rather than discussed under the separate 
phases. 


Materials, Equipment, and Specimens 


Materials 


Cement 


Type If cement was used throughout the investigation and was obtained 
from the same mill but in three different shipments. The chemical and physi- 


cal data for each shipment, designated by the numbers RC-183, RC-186, and 
RC-233, are shown below: 


Stage 1 Stage 2 Stage 3 


Chemical Properties 4 
Constituents, 

22.3 22.5 22.2 
A1903 4.3 5.3 4.8 
Fe 203 3.5 3.2 3.7 
CaO 63.2 64.0 63.9 
MgO 3.3 3.0 3.1 
S03 1.5 1.5 1.5 
Ignition loss 0.55 0.68 0.52 
Insoluble residue 0.24 0.15 0.19 
Na2O 0.17 0.21 0.13 
K90 0.48 0.34 0.41 


Total alkalies as Na20 


0.49 


PRESSURE GROUTING 
Calculated compounds, 
C35 45 46 50 
C25 2 30 26 
C3A 9 6 
C4AF 12 10 ll 
CaSO04 3 
Physical Properties a 
. Normal consistency 25.2 24.2 25.0 
Setting time, Gilmore, hr-min: 
4 Initial 4:00 3325 3315 
Final 6:20 6325 4:15 
“| Autoclave expansion, » 0.09 0.10 
Air content of mortar, % 4.9 7.3 3.5 
sompressive strength of mortar, psi: 
3 days 260 7 
Fineness, cm 
Wagner 1805 720 1840 
Blaine 3170 2965 3080 
Sieve analysis, dry: 
No. 50 (x% retained) 0.2 0.0 0.0 
Noe LUO retained) 0.1 0.1 0.0 
No. 200 (% retained) 264 4.3 1.3 
Fly Ash 
The fly ash used in all three stages came from the Chicago region but in 
different shipments. A chemical analysis was not made of the fly ash used in 
stages 1 and 2. The chemical analysis of the fly ash used in stage 3 is shown 
in the following paragraph. The physical characteristics of the fly ash used 
in the respective stages are listed in the following tabulation: 
Opening ____Per Cert Retained 
Sieve Nos — = Stage 1 Stage 2 Stage 3 
50 0.0117 0.4 0.2 0.9 
100 0.0059 1.3 lel 0.9 
200 0.0029 6.2 6.0 245 
Fineness, 
Blaine 2795 2795 —_ 
Fisher -- -- 3335 
Specific gravity 
Other Mineral Admixtures 
The following mineral admixtures were used in stage 3 only: 
a. Slag, water-quenched, ground, laboratory No, RC-216, 
b. Pumicite, laboratory No. AD-6. 
c. Opaline shale, calcined, laboratory No. AD-13. 
Data on these materials and on the fly ash used in stage 3 follow: 
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Opaline 
Fly Ash Slag Pumicite Shale 
AD-3 RC-216 AD-6 AD-13 
Chemical Properties . 
Constituents, % 
Moisture loss 0.20 oe 0.83 1.2 
S$i02 47.2 35.8 68.8 70.1 
A1203 19.5 15.8 14.8 19.3 
Fe203 18.2 1.2 1.4 
Mn203 0.07 0.62 0.05 0.05 
P205 0.23 a= 0.02 0.14 
CaO 5.3 35.5 0.65 0.2 
MgO 1le2 10.4 0.33 0.35 
Sulfide sulfur 0.05 -- 0.00 0.00 
S03 262 0.1 0.03 0.10 
Ignition loss 0.8 0.23 3.9 1.26 
Insoluble residue 70.4 0.52 94.6 86.9 
Na20 (gravimetric) 0.69 -- -- -- 
K20 (gravimetric) 1.16 -- -- -- 
Total Na20 (gravimetric) 1.45 -- -- -- 
Na20 (flame) 1.62 0.11 1.38 0.40 
K20 (flame ) 1.98 0.65 4.96 0.61 
Total Nag (flame) 2.92 0.54 4.64 0.80 
CHC13 sol -- 
Total carbon 0.43 0.54 0.02 0.07 
Opaline 
Fly Ash Slag Pumicite Shale 


AD-3 RC-216 AD-6 AD-13 


Physical Properties 
Fineness, cm2/g: 


Wagner -- 3,430 -- -- 
Blaine -- 5,795 -- -- 
Fisher 3,335 5,470 4,640 15,020 
Specific gravity 2.48 2.85 2.35 2.35 
Sieve No., and sieve opening (in.) 
50 0.0117 0.9 0.8 0.0 0.1 
100 0.0059 0.9 0.4 0.4 0.9 
200 0.0029 225 1.3 1.5 13.9 


Chemical Admixtures 


The intrusion aid causes slight expansion of the grout and is purported to 
prevent agglomeration of the solids and aid penetration of fissures. It was 
used in all tests at the rate of one per cent of the weight of cement or cement j 
plus fly ash whenever fly ash was used in the grout. 7 
Calcium lignosulfonate, cls, was added for the same purposes as intrusion 
aid. It did not, however, cause expansion of the grout. It was used in all 
cases at the rate of 0.23 per cent of the weight of cement. 


a 


» 


PRESSURE GROUTING 
Grouting Equipment 


The grout pump was a single-cylinder, reciprocating 3-3/4- by 2-1/2- by 
5-in, steam pump with rubber piston and valves, air-driven for these tests. 
The grout mixer was a paddle-type machine powered by a motor from an 
electric drill. A photograph of the grouting assembly, including a typical 
specimen, as used in stage 1 of the program is shown on Fig. 1. An important 
change in the equipment setup was made at the conclusion of stage 1 and was 
used throughout the remainder of the program. This change is shown sche- 
matically on Fig. 2 and consisted of connecting the delivery line and a return 
line with a corresponding control valve immediately adjacent to gage A (gage 
O on Fig. 2). It was suspected that the velocity of flow in the riser of the de- 
livery line would be so low in some cases as to allow settlement of the 
coarser particles and thus result in a variable grout passing through the 
fissure. This supposition was confirmed by observing the flow of grout 
through a vertical lucite tube. This difficulty was largely overcome by open- 
ing the valve on the return line, as shown on Fig. 2, sufficiently to insure a 
relatively constant rate of grout flow through the riser under all grouting con- 
ditions. In addition, three additional gages were installed along the specimen 


so that variations in pressure gradient could be determined during grouting 
operations, 


GAGE A 
TEST SPECIMEN 


INPUT LINE 


PRESSURE CONTROL VALVE ~~ 
AIR LINE” 
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Specimens 


A specimen consisted of a pair of matching prisms of high-strength rein- 
forced concrete. Each half was 7 in. wide by 3 in. thick by 52 in. long. Two 
prisms were cast at the same time in the same form with a piece of 1/4-in. 
plate glass as a separator between them to assure matching surfaces. Half 
of the prisms were without pipe connections. Half of those in the first stage 
tests contained two 3/4-in. pipe nipples 48 in. apart embedded along the 
center line normal to and in contact with the glass separator. Half of those 
in the second and third stage tests contained five nipples 12 in. apart. 

A pair of matching prisms was separated by a ribbon of shimstock 1/2 in. 
wide of the required thickness around the edges and ends between the halves 
to form a fissure 6 in. wide and 48 in. long between the end nipples. The two 
halves were held together by C-clamps. The specimens were pumped ina 
horizontal position for the following reasons: 


a. The most dangerous conditions in grouting under a dam are caused by 
extensive open passages such as horizontal bedding planes. It was de- 
sired to learn something about grout flow and pressure distribution 
under such conditions. 

b. It was hoped to learn something concerning the validity of the as- 
sumption that excess mixing water can be forced out of the grout by the 
application of pressure. It was believed that if this did occur a good 
bond would be made with the top as well as the bottom half of the speci- 
men; if not, bleed water would prevent bond with the top half. 
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Grout Pumping Tests 


General Procedures 


Treatment of Grout 


Dry materials were all sieved through a 30-mesh sieve and the mixed 
grout was passed through a 30-, 50-, or 100-mesh sieve depending upon the 
water-cement ratio. The thinnest grouts were passed through a 100-mesh 


sieve, thicker grouts through a 50-mesh sieve, and the thickest grouts through 
a 30-mesh sieve. 


Grout Injection 


Grout injection was started by pumping water through the specimen 
followed by thin grout which was gradually thickened by the addition of solids. 
The desired pressure was maintained on the specimen by operating the bypass 
valve so that only a controlled portion of the total flow of grout was shunted 
through the specimen. The pump was operated at a speed of approximately 
72 strokes per minute for the major portion of the tests. At the lowest water- 
cement ratio for a given test more than one hour of pumping was usually re- 
quired to discharge one cubic foot of grout through the fissure of a specimen. 
Grout flow stopped at any given pressure when the grout became too thick. 

In stage 1 operations, pressure was maintained on specimens for 10 to 15 

min after all flow had stopped with the valve at B (Fig. 1) open at all times. 
Water continued to drip very slowly from the discharge line during final ap- 
plication of pressure, being squeezed out after separating from the grout film 
in the specimen. When gage B assembly was removed the short pipe nipple 
cast in the test specimen was invariably found to be full of hardened grout into 
which a pencil could be forced only with difficulty. This was the case regard- 
less of the consistency of the grout being pumped. In stage 2 operations the 
valves at the entrance and discharge ends of the specimen were closed as 
soon as all flow had stopped and the pressure within the specimen was allowed 
to equalize insofar as internal conditions would permit. Pumping for stage 3 
was in most instances conducted as in stage 1 although for some tests the 


valves were closed at cessation of pumping on both ends of the specimens, as 
in stage 2. 


Test Conditions 


The study was conducted in three separate stages as follows: 


a. The first stage furnished data on the lowest water-cement ratios of 
grouts that could be pumped through fissures of 0.01-, 0.02-, and 0.03- 
in. thickness at 100 psi pressure, with the following grout conditions: 
(1) neat cement; (2) cement plus fly ash; (3) cement plus fly ash plus 
the admixture, intrusion aid. A study of consistency, bleeding charac- 
teristics, and setting times of the various grouts was also made. 

b. The second stage provided data on grout penetration obtained at pump- 
ing pressures of 25 and 50 psi using the three grout conditions of the 
first stage and the additional condition of cement plus intrusion aid. 
Tests were also made of neat-cement grout plus calcium lignosulfonate, 
hereinafter abbreviated as cls, at 25, 50, and 100 psi, and of neat- 
cement grout plus intrusion aid at 100 psi. 
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c. The third stage consisted of: pumping tests through an 0.03-in. thick 
crack at 50 psi using the seven grout conditions listed below; tests of 
the grouts themselves for consistency, bleeding, setting time; and tests 


to determine the apparent solubility of hardened grout films in distilled 
water. 


Grout Conditions 


(1) Neat cement (5) Cement plus slag 
(2) Cement plus fly ash (6) Cement plus pumicite . 
(3) Cement plus cls (7) Cement plus opaline shale 


(4) Cement plus fly ash plus cls 


Some of the above materials were used in various proportions of one 
to the other so that a total of 17 different grout mixtures were pumped. 
These are described in detail later in the report. 


Observations were made during and after pumping for: 


a. The effect of specimen surface on the water-cement ratio and pump- 
ability of the grout (stage 1 only). 

b. Lowest water-cement ratio of grout that would penetrate the three 
different fissure thicknesses. 

c. Pressure drop along the specimen. 

d. Rate of grout flow at a given pressure, water-cement ratio, and crack 
thickness. 

e. Relationship between pumping pressure and lowest water-cement ratio 
grout that a fissure would take. 

f. Relationship of particle size of grout solids to the width of crack that 
could be grouted. 

g. Effect of chemicals (intrusion aid and cls) on the penetration qualities 
of grouts. 

h. Quality of hardened grout films, from visual observation. 

i. Strength of grout film as judged by shear tests (for stage 2 only). 

j. Special observations and tests were made in stage 3 for consistency of 
the grout at each change in water-cement ratio by measuring change in 


unit weight of the grout and the torque imparted to a piano-wire con- 
sistency meter. 


Discussion of Pumping Tests 


Ninety-seven pumping tests were made. The detailed data, because of 
their voluminous nature, are not tabulated in this report. They are presented 
and analyzed in the figures and discussion. In most instances the lowest 
water-cement ratio shown for a given grout mixture and test condition was 
the lowest water-cement ratio that could be pumped. It is considered, how- 
ever, that the water-cement ratio at which the flow equals 0.1 cu ft per min 
(6 cu ft per hr) is the stiffest grout practical for field pumping. 

Throughout this report water-cement ratio is expressed by weight of water 
to cement or to cement plus other solids when materials such as fly ash were 
added to the grout. All references in text or figures to cement-fly ash grout, 
or cement-slag grout, etc., refer to grout having equal parts of cement to fly 
ash or other material by weight unless specifically stated otherwise. 


. 
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Effect of Surface Texture of Specimen on Grout Penetration 


Pumping tests 2 through 12 were made using specimen surfaces from 
which the glaze, due to casting against the glass, was removed by light rub- 
bing with a carborundum stone. Tests 13 through the end of the program were 
made using specimen surfaces from which the glaze had not been removed. 
The effects of type of specimen surface, fineness (maximum grain size) of the 
grout, and crack thickness on the water-cement ratio of the grout are shown 
on Fig. 3 and discussed below: 


a. The combination of smooth surfaces and grout strained through a No, 50 
sieve permitted the use of lower water-cement ratios for the 0.03- and 
0.02-in. fissures than for the combination of rubbed surfaces and grout 
strained through a No. 30 sieve, for both neat grout and grout contain- 
ing fly ash. 

The combination of smooth surfaces and grout strained through a No. 
50 sieve permitted the passing of neat grout with a water-cement ratio 
of 1.34 by weight through the 0.01-in. opening as contrasted with a 
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water-cement ratio of 2.06 with the rubbed surfaces and neat grout 
strained through a No. 30 sieve. 

c. Grout containing fly ash could not be pumped through the 0.01-in. open- 
ing regardless of the type of surface or when strained through the No. 
50 sieve even with a water-cement ratio of 4.30, the highest water- 
cement ratio used in the tests. The reason for the poor pumpability of 
the fly ash grout is discussed in the following paragraph. 

d. The use of one per cent of intrusion aid by weight of the cement plus fly 
ash permitted a slight reduction in water-cement ratio when grouting 
the 0.03-in. fissures with rubbed surfaces as compared to the grout 
containing fly ash and no intrusion aid, but did not help when grouting 
the 0.02- and 0.01-in. fissures with rubbed surfaces. 

e. The use of intrusion aid permitted a slight reduction in water-cement 
ratio when grouting the smooth-surfaced 0.02- and 0.03-in. fissures but 
was ineffective in promoting penetration of the 0.01-in, fissures. 


The physical data on the materials used in stage 1 showed the cement to 
have a greater specific surface (Blaine 3172 cm2 per g, 97.3 per cent passing 
the No. 200 sieve) than the fly ash (Blaine 2795 cm2 per g, 92.1 per cent pass- 
ing the No. 200 sieve), and showed that the fly ash contained somewhat more 
than 0.4 per cent of material larger than 0.0117 in. (No. 50 sieve). This grain 
size exceeded the 0.01-in. opening of the smallest fissure and caused a grid 
of oversize grains to build up quickly during pumping which blocked the 
passage of the grout through the fissure regardless of whether or not the 
grout contained intrusion aid. 


Influence of Crack Thickness on Grout Penetration 


Fig. 4 is a bar graph showing the water-cement ratio for each type of grout 
at the stage where it just failed to penetrate the three fissure widths at the 
pumping pressures used. Grouts of only slightly higher water-cement ratios 
than those indicated penetrated the fissures. For neat-cement grout pumped 
at 25 psi the lowest water-cement ratio that would penetrate an 0.01-in. crack 
exceeded 2.67, when this pressure was increased to 50 psi the water-cement 
ratio exceeded 1.33 but an increase to 100 psi did not permit further re- 
duction in water-cement ratio. 

Little practical difference was noted between the water-cement ratios of 
the neat grouts that would penetrate the 0.02- and 0.03-in. cracks. Pumping 
pressure, within the range investigated, also had little effect on the ability of 
grouts of various water-cement ratios to penetrate the 0.02- and 0.03-in. 
cracks. 

Grouting of any of the 0.01-in. cracks with grouts containing fly ash was 
found impracticable for the reasons previously stated. 

The use of fly ash was of no benefit in grouting the 0.02-in. crack, since it 
permitted no practical reduction in water-cement ratio over neat grout. 
Pumping pressure, within the range used, was of little value in improving the 
penetration of the 0.02-in. crack with a cement and fly ash grout. 

Grouts of cement and fly ash that would penetrate the 0.03-in. crack were 
of slightly lower water-cement ratio than the neat grouts that would penetrate 
cracks of similar width at 25 and 50 psi, but practically the same at 100 psi. 

Use of intrusion aid in the neat-cement grout permitted a slight reduction 
in water-cement ratio of the grouts that would penetrate both the 0.02- and 
0.03-in. cracks at the pressures tested. 
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Cls was used only at 50-psi pressure on the 0.02-in. crack with neat 
cement and with cement plus fly ash but appears to be as effective as in- 
trusion aid with respect to grout penetration. 

In the above discussion the effect of maximum grain size should be kept in 
mind as well as the fact that Fig. 4 indicates water-cement ratios at refusal, 


whereas a rate of flow of 6 cu ft per hr is the stiffest grout practical for field 
pumping. 


Pressure Drop Along the Specimen 


The pressure gradient for the 0.02- and 0.03-in. cracks during pumping 
was approximately a straight line with zero pressure at the exit end of the 
specimen (gage 4), maximum pressure at the entrance end of the specimen 
(gage 0), and intermediate pressures at gages 1, 2, and 3. The static 
pressure, after pumping was stopped when the grout was thick (water-cement 
ratio of less than 0.5) and valves at both entrance and exit ends were closed, 
did not always equalize from end to end. This may have been due to a thixo- 
tropic condition developing in the thick grout immediately upon cessation of 
movement, and causing a blocking of the specimen between gages. For the 
0.01-in. crack, the pressure drop during pumping did not follow a straight 
line. Fig. 5 is a series of curves of pressure gradients that developed when 
neat grouts were pumped at various water-cement ratios through an 0.01-in. 
crack at 25 and 50 psi. It will be noted that the gradient was fairly uniform 
along the whole specimen length while water was being circulated through it; 
however, as the grout was thickened flow diminished and the pressure gra- 
dient between the first two gages (first 12 in. of length) increased. This would 
apparently indicate that in the very fine seams as the grout thickened and the 
flow diminished the hydraulic pressure fell off rapidly with distance from the 
intrusion point. Conditions of pressure gradient and static pressure after 
pumping are shown on Figs. 6-14 for crack thicknesses of 0.02 and 0.03 in, 

These figures show that in cracks with smooth walls and thicknesses of 
0.02 in. or more the conditions of flow and pressure are similar to those of 
water. Irregularities in the curves are due to faulty operation of the gages oc- 
casioned by plugging. Water-cement ratio values at which the grouts would 


act like water, relative to pressure distribution, might be different in nature 
where fissure walls are rough. 


Rate of Flow at Given Pressure, Water-Cement Ratio and Crack Thickness 


Fig. 15 is a bar graph showing the amount (in cubic feet per hour) of vari- 
ous grouts, all with a water-cement ratio of 0.5, that could be pumped through 
0.02- and 0.03-in. cracks at the three pumping pressures employed. Fig. 16 ° 
presents the same data in curves, omitting the data on the cement fly ash 
plus intrusion aid grout, which are anomalous and inconclusive. 

Increasing pressure caused proportionately increasing flow. Apparently 
the flow of the cement fly ash grout was greater at equivalent pressures and 
crack openings than that of neat grout, 

The use of intrusion aid in the neat grout increased flow at equivalent 
pressures and crack thicknesses over plain neat grout and caused greater 
flow than did the use of fly ash. Data on the effect of intrusion aid on the 
cement fly ash grout were inconclusive, 

The rate of flow provides a realistic measure of the viscosity of a grout. 
Figs. 15 and 16 show that the rate at which a grout can be injected into a 
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crack is influenced by pressure and composition, Rate of flow is an index to 
the distance a grout can be pushed at a given pressure before friction causes 
the flow to cease. Under conditions dictating the use of low-pressure grout- 


ing, greater density of drilling than ordinary or the use of fluidifiers is indi- 
cated. 


Relationship Between Pumping Pressure, Water-Cement Ratio, and Crack 
Thickness 


An examination of Fig. 4 shows that the crack thickness had a great effect 
on the lowest water-cement ratio grout that could be pumped through a crack, 


The pressure had little practical effect, since water-cement ratios lower than 
0.5 are seldom employed. 


Relationship Between Particle Size of Grout Solids and Width of Crack 
Penetrated 


It will be noted that the lowest water-cement ratio neat grout that could be 
forced through the 0.01 crack (Fig. 4) was approximately 1.33 by weight. This 
is a dilute material of dubious quality and is believed to be neither durable 
nor impermeable, Grout containing fly ash could not be forced through the 
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0.01-in. crack even at considerably higher water-cement ratios. Examination 
of the physical data on the cement and fly ash shows that the maximum grain 
size of the cement was between 0.01 and 0.006 in., corresponding to the 50- 
and 100-mesh sieves with 0.1 per cent retained in this size range. The maxi- 
mum grain size of the fly ash was somewhat in excess of 0.01 in. These re- 
sults indicate that stoppage of the crack occurs when the crack-opening: 
grain-size ratio is 1.7 even with quite dilute grouts, and that with fly ash, 
where the maximum grain size exceeds the crack opening, grouting is im- 
possible regardless of the water-cement ratio. The critical crack-opening: 
grain-size ratio probably varies with different materials, and it is believed 
that a safe value should exceed 3.0. 

Alfred Machis(1) found in the grouting of sand that the ratio of pore diame- 
ter to grain size necessary to permit penetration of cement slurry was at 
least five. A. F. Taggart, (2) in his handbook, shows a graph that indicates 
the ratio of filter opening to maximum particle diameter for filter-cake for- 
mation. The fly ash had a maximum grain diameter somewhat larger than 
250 microns, From the curve cited above, a filter-opening: grain-size ratio 
somewhat greater than three is indicated as necessary to prevent formation 
of a filter cake. The shape of the opening being grouted should influence the 
ratio of particle diameter to thickness. A slot-shaped opening, as used in the 
tests being reported, because of its relatively infinite width and no restricting 
sides, should permit a smaller ratio of particle-size diameter to opening to 
pass grout than a square or rounded opening of the same thickness. 


Effect of Intrusion Aid and cls on Penetration Characteristics of Grouts 


Fig. 4 shows that when intrusion aid was added to neat grout when grouting 
an 0.02-in. crack the resulting water-cement ratio was reduced 0.03 at 25- 
psi, 0.08 at 50-psi, and 0.17 at 100-psi pumping pressure, or an average re- 
duction of 0.09 for all three pressures. When intrusion aid was used in 1 
cement: 1 fly ash grout, no reduction in water-cement ratio but rather an in- 
crease of 0.43 resulted in grouting the 0.02-in. crack at 25 psi. However, an 
0.02 reduction in water-cement ratio resulted at 50 psi and an 0.09 reduction 
at 100 psi. 

The use of intrusion aid in grouting the 0.03-in. crack with neat grout al- 
lowed a reduction in water-cement ratio of 0.07 at 25 psi, 0.11 at 50 psi, and 
an unknown amount at 100 psi since the cement plus intrusion aid grout at 
this pressure was not thickened to refusal. The use of intrusion aid in grout- 
ing the 0.03-in. crack with cement fly ash grout apparently caused an in- 
creased water-cement ratio of 0.03 at 25 psi, a decrease of 0.01 at 50 psi, 
and a decrease of 0.03 at 100 psi. 

Fig. 4 also shows that the use of cls in neat grout permitted a reduction of 
0.01 at 25 psi, of 0.12 at 50 psi, and of 0.09 at 100 psi, or an average re- 
duction of 0.07 in the minimum water-cement ratio that would penetrate an 
0.02-in. crack. In grouting the 0.03-in. crack, the average reduction in water- 
cement ratio was 0.04 for the three pumping pressures used. The use of cls 
in cement fly ash grout apparently allowed a reduction of 0.07 at 25 psi, of 
0.03 at 50 psi, and an increase of 0.01 in the water-cement ratio at 100 psi 
when grouting the 0.02-in. crack, In grouting the 0.03-in. crack with cement 
fly ash grout a slightly increased water cement ratio was required for pene- 
tration at all three pressures used. 

For penetration of a given crack thickness, intrusion aid seemed to be a 
little more efficient in lowering the water content of grout than cls. That is, 


| 
7 


ASCE PRESSURE GROUTING 1731-19 


intrusion aid was slightly more efficient than cls as a fluidifier. However, 
little appears to be gained from use of either material except where it might 
be desirable to inject low water-cement ratio grout under low pressure. 


Quality of Hardened Grout Films from Visual Observation 


The two halves of the grouted specimen were opened and examined approxi- 
mately 24 hr after pumping. The following tabulation summarizes these 
observations with the evaluations arranged in order of increasing water- 
cement ratios for each type of grout. 


Bond Judged 
Water-cement from 
—___—Type Grout _Ratio, wt _Appearance  _Appearance _ 
Neat 0.32-0.38 Good Good to poor 
0.46 Good Fair 
0.45 Fair Poor 
0.51 Fair Poor 
1.33 Poor, stringy Poor 
1.34 Poor, stringy Poor 


Neat + intrusion aid 0.34 Good Good 
0.34 Good Fair 
0.35 Good Good 
0.36 Good Poor 
0.40 Good Fair 
0.43 Good Fair 
4.0 Poor, stringy Poor 


Neat * cls 0.36 Good Fair to poor 


1:1 (cement:fly ash) 0.36 Good Poor 
0.38 Good Poor 
0.40 Good Poor 
0.45 Fair, stringy Poor 
0.49 Good Poor 
4.0 Poor, stringy Poor 
4.3 Poor, stringy Poor 


1:1 * intrusion aid (cement: 0.33 Good Poor 
fly ash * aid) 0.37 Good Poor 
0.40 Good Poor 
0.43 Good Fair 
0.43 Fair Poor 
3.0 Poor, climped Poor 
4.3 Poor, stringy Poor 


1:1 + cls (cement:fly ash + 0.035 to Good Good 
els) 0.038 


1.5:1 (cement:fly ash) 5 Fair, stringy Poor 
2:1 (cement:fly ash) : Fair, stringy Fair 
1.5:1 (cement:fly ash + cls) Good Good 
2:1 (cement:fly ash * cls) Good 
1:1 (cement:slag) Good 


1.5:1 (cement: slag) 3 Good 
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Bond Judged 
Water-cement from 
Type Grout Ratio, wt Appearance Appearance 
2:1 tement:slag) 0.38 Good Poor 


1:1 (cement: pumi cite) 0.56 Fair (soft) Fair 
1.5:1 (cement: pumicite) 0.50 Good Fair 


2:1 (cement:pumicite) 0.50 Good Poor 


1:1 (cement:opaline shale) 0.60 Fair (soft) Poor 
1.5:1 (cement:opaline shale) 0.55 Good Poor 


2:1 (cement:opaline shale) 0.55 Good Poor 


A grout film designated “good” in the tabulation was one which was hard to 
the fingernail, filled the cavity completely, and showed “none” to a “moderate” 
amount of bleeding as evidenced by channels cut by “bleed water” as it moved 
under slight residual pressure from the pumping process, or under gravity, 
toward points of lesser pressure or lower elevation. Fig. 17 shows a film 
designated “good”. A film was designated “fair” when the grout could be 
scratched with the thumbnail or was traversed by numerous bleeding channels. 
It was called “poor” when it was quite soft or when the film occurred in 
lenses or streaks and could not be considered a continuous sheet of material. 
Fig. 18 shows a close-up of a poor-quality film. 

Bond was considered “good” when approximately 40 per cent or more of 
the film adhered to the upper slab when the specimen was opened; “fair” when 
the adherence was approximately 25 to 40 per cent; and “poor” when the ad- 
herence was less than approximately 25 per cent. 

It will be noted that the appearance of all the grouts having a water-cement 
ratio of 0.6 or less was “good” to “fair,” and that where “poor” grout was en- 
countered the ratio was 3.0 or higher. There are no intermediate values with 
the water-cement ratios ranging from 0.6 to 3.0, but it is only logical to as- 
sume that the quality of such films would have been “fair” to “poor.” A di- 
rect comparison of the quality of different types of grout having the same 
water-cement ratio is not possible from the present data. However, judging 
from the bond evaluations, the use of intrusion aid appeared to be of some 
slight benefit in increasing adherence of the film to the top slab; but for grout- 
ing a horizontal opening between two nonabsorbent surfaces where a high- 
quality grout would be required, which would assuredly adhere to the top as 
well as the bottom of the cavity, a water-cement ratio not exceeding approxi- 
mately 0.46 and neat grout should be used. The use of fly ash, judging from 
these qualitative data, may have had an adverse effect on the bond between 
top and bottom slab, although the quality of the films, from visual obser- 
vation, was usually good. There was an apparent tendency toward clumping 
of the solids in some cases. 

Cls in the cement fly ash grouts tended to make them softer at 24 hr, but 
appeared to counteract to a large extent the tendency to agglomerate, and 
improved the bond to the top slab. 

The use of slag, pumicite, and opaline shale tended to improve the ap- 
pearance of the grout films to a marked extent. Visual evidence of bleeding 
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was entirely lacking with all the pumicite and opaline shale grouts and with 
the 1 slag: 1 cement grout. The 1.5 cement: 1 slag grout showed slight bleed- 
ing and the 2 cement: 1 slag grout showed moderate bleeding. All three slag- 
grout films were hard to the fingernail at 24 hr; however, the 1 cement: 1 
pumicite grout could be dented by a fingernail at 24 hr and the 1 cement: 1 
shale grout could be slightly scratched by a fingernail at 24 hr. Bonding to 


the top slab did not seem to be improved by use of slag, pumicite, or opaline 
shale, 


Strength of Grout Films Judged by Shear Tests 


An attempt was made, during stage 2 of the tests, to measure the bond 
strength developed between top slab and the grout by making shear tests at 
120 days age on sections of grouted specimens cut to 10-in. lengths. Three 
0.01-in. cracks were grouted at 50 psi with various grouts as listed below, all 
with a water-cement ratio of 0.45 and with pressure maintained on the speci- 
mens for 15 min after pumping: 


Type Grout Test No. Remarks 
Neat + intrusion aid 46A 
Neat 48 
1 cement: 1 fly ash 49 Came apart in handling 
1 cement: 1 fly ash 54 Came apart in handling 


1 cement: 1 fly ash + aid Came apart in handling 


The following results were obtained: 


Specimen Section Shear Resistance, psi 
46(1) 50 
46(2) Came apart in handling 
46(3) 180 
46(4) 70 
48(1) 65 
48(2) 30 
48(3) 50 
48(4) 60 


Test values were quite low for the neat-grout and neat-grout plus aid speci- 
mens tested, indicating very little bond. Those specimens containing fly ash 
did not develop enough bond to permit making a test. 


Tests of Grout Characteristics 
Consistency 


A satisfactory device for measurement of consistency is highly desirable 
and would enable the grouting crew to suit the quality and “pumpability” of 
the grout to the grouting conditions encountered. Consistency measurements 
were made on the grout itself by means of a piano-wire, pendulum-type 
viscosimeter, Fig. 19, by a flow cone, and by a unit-weight device. 

The viscosimeter was developed by Professor R. E. Davis of the University 
of California, Berkeley, and the Bureau of Reclamation. The viscosity or 
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Pig. 19. Pilanowire viscosimeter 


consistency of a slurry is measured by placing a sample in a shallow metal 

container mounted on a motor-driven turntable. A wire spider is suspended 
from a piano wire and is submerged in the sample of grout. A weight holds 

the piano wire taut. Torque is imparted to the spider and piano wire by the 

turning container of grout. The number of degrees of torque varies with the 
consistency of the sample being tested. 

The flow cone is simply, as its name implies, a cone with a specified ori- 
fice from which a known volume of liquid flows in a certain length of time, 
depending upon its consistency. The flow cone contained 1725 ml which were 
discharged through a pipe orifice, 1-1/2 in. long by 1/2-in. inside diameter, 
in the bottom of the cone. Flow was timed with a stop watch. In some of the 
tests the discharge pipe was modified to 6 in. in length and 3/8-in. inside 
diameter. 

The unit-weight device for consistency determinations was designed for 
use with neat-cement grout and consisted of a measure and specially graduat- 
ed beam for translating unit weight into water-cement ratio. 

The piano-wire torque meter provided the most satisfactory means of 
measuring consistency in the laboratory. It was calibrated against the 
Stormer viscosimeter in connection with another investigation so that the 
readings obtained could be translated into standard units of viscosity (poises). 
This calibration is shown on Fig. 20. Consistency meter readings at fixed 


water-cement ratios for grout combinations used in the third-stage work are 
tabulated below: 
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Neat cement 

1 cement:1 fly ash 
2 cement:1 fly ash 
Cement * cls 
cement:l1 fly ash + cls 
cement:1 fly ash * cls 
cement:1 slag 

cement:1 slag 

cement:1 pumicite 
cement:1 pumicite 
cement:1 opaline shale 
cement:1 opaline shale 


the solids in the grouts. 


52T insofar as practicable. 


A water-cement ratio of 0.4 was used in all cases except for the grouts 
containing opaline shale, which were too thick for use at 0.4, so 0.6 was used 
for these. Placing the grouts on approximately the same water-cement ratio 
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Consistency, deg Torque _ 


__Water-cement Ratio, wt 


350 
332 
337 
374 
258 
263 
378 
355 


* Exceeded the capacity of meter. 
#* Not fluid grout, too stiff to test. 


Bleeding 


026 


For field work the torque meter might be too delicate since it would require 
a special indoor working area free of vibration. 
For field work specially constructed flow cones with volume and discharge 
orifice integrated to produce sizable differences in readings with small chang- 
es in viscosity, or unit weight measures with the unit weight of the grout ex- 
pressed in grams per gallon would probably be satisfactory. 
Whatever the method used for measuring viscosity, it should correlate 
with pumpability as measured in cubic feet per hour in order to have any 
practical meaning. Examination of Fig. 21 shows good correlation between 
torque and pumpability (flow through the 0.03-in, crack at relatively constant 
pressure, approximately 50 psi) for several combinations of materials. 
There appeared to be little correlation between torque and flow when the 
grout combinations were considered collectively. This apparent lack of simi- 
larity of pumpability for grouts of equal viscosity through a thin crack is 
probably due to differences in grain-size distribution and particle shape of 


Two kinds of tests were carried out to determine the extent of bleeding of 
the grouts. The first consisted simply of pouring 500 ml of grout into a 1000- 
ml graduate and observing the separation of water with passage of time. The 
water-cement ratio was varied and measurements were made on neat grouts, K 
1 cement: 1 fly ash grouts, and 1 cement: 1 fly ash grouts with 1 per cent in- 
trusion aid. Fig. 22 shows data obtained after one hour. Most of the bleeding 
had occurred at the end of one hour; no additional bleeding occurred after two 
hours. Use of intrusion aid in cement fly ash grout appeared to reduce slight- 2 
ly its bleeding below that of neat and cement-fly ash grouts, 

The second kind of bleeding test was performed according to ASTM C 243- 
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basis regardless of solid material composition permitted evaluation of the ef- 
fects of the constituents, other than water, on their bleeding characteristics. 


The combinations used and results obtained are shown in the following tabu- 
lation: 


wW/C, 0.4, wt W/C, 9.6, wt 
Rate Rate Capacit 
em3/em2/ em3/cm2/ 


Neat cement 80 37 

1 cement: y ash 65 38 

2 cement: ash 45 

Cement. + c f 6 

cement: ash * cls 20 

cement: y ash * cls 36 

cement:1 s 23 

cenent:1 4 33 

cement:1 pumicite 5 18 

cement:1 pumicite 14 -- 
cement:1 opaline shale * 34 5 
cement:1 opaline shale 72 18 


Note: The rate indicates how rapidly the cement grains settle and the 
water rises and collects on the top surface of the paste. The 
capacity indicates the amount of water, per unit volume of paste, 
that collects on the surface during the test. 


* Not fluid, too stiff to test. 


Cls appeared to decrease sharply the rate and amount of bleeding when 


used in neat grout. It appeared to act in a similar manner in the grouts with 
fly ash but to a lesser extent. Fly ash itself appeared to slow down the rate 
of bleeding but did not affect the total amount of bleeding. Slag appeared to 
lessen both the rate and amount of bleeding to a small extent. Both the pumi- 
cite and opaline shale sharply decreased the amount of bleeding and, to a 
somewhat lesser degree, the rate of bleeding. 


Setting Time 


Setting times of grouts containing the same constituents, except for the ad- 
dition of intrusion aid in two cases, as used for the second series of bleeding 
tests, and made with several water-cement ratios, were determined. The de- 
terminations were made by use of the 1-mm Vicat needle and smal] samples 
of grout placed in shallow wide-mouthed vials. The data obtained are shown 
in the tabulation following. 

Setting time is obviously important when grout is expected to cut off perco- 
lating water. In grouting by stages, it also controls the time that must elapse 
before holes can be cleared out and drilling resumed. 

Setting times could not be determined exactly because of hours-of-work 
limitations. Cls appeared to lengthen setting time. Setting time was also 
influenced by water content. It is conceivable that grouts could have water 
contents so high that the cement grains would be separated to such an extent 
that setting and development of strength could never occur. 


a. 
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Two 50-ml portions of grouts, listed in Fig. 23, were placed in stoppered 
200-ml flasks and one was cured for 7 days, the other for 90 days. At the end 
of the curing age, the samples were leached with 50 ml of distilled water for 
periods of 2, 7, 28, and 90 days. Leach water was then tested for pH, total 
solids, and calcium ion. The amount of total solids and calcium removed by 
leaching is shown graphically in Fig. 23. 

Unfortunately reaction occurred between the Pyrex glass flasks in which 
the grouts were cured and the alkalies in the grouts, causing many to crack 
after several weeks storage permitting the water to leak out. No test results 
are plotted from specimens which were in the cracked flasks. 

Conclusions drawn from the unspoiled tests are: 


Water-cement Approximate 

Neat cement 0.33 7-1/4 
Neat cement 0.4 Exceed 4, under 18 
Neat cement 0.8 Exceed 4, under 19 
Neat cement 2.06 233 

fi 1 cement:1 fly ash 0.4 Exceed 6, under 22 

: 1 cement:1 fly ash 0.44 16-1/2 
1 cement:1 fly ash 0.80 Exceed 5, under 21 
1 cement:1 fly ash + 1% intrusion aid 0.36 20 
1 cement:1 fly ash + 1% intrusion aid 4.32 Not set in 336 
Cement + cls 0.4 Exceed 7, under 23 
Cement + cls 0.8 72 
1 cement:1 fly ash * cls 0.4 Exceed 22, under 46 
1 cement:1 fly ash + cls 0.8 Exceed 22, under 46 
1 cement:1 slag 0.4 Exceed 6, under 22 
1 cement:1 slag 0.8 30 
1 cement:1 pumicite 0.4 Exceed 4, under 20 
1 cement:1 pumicite 0.8 Exceed 6, under 22 
1 cement:1 opaline shale 0.6 Exceed 3-1/2,under 18 
1 cement:1 opaline shale 0.8 30 


Solubility 


Grouts exposed to leaching after only 7 days curing appeared to be 
somewhat more susceptible to leaching than those cured 90 days. 
About one-half of the dissolved material that leached out of the 0.4 
water-solids ratio grouts was calcium, whereas about one-third of the 
dissolved material from the 0.8 water-solids ratio grouts was calcium. 
Increasing the water-solids ratio from 0.4 to 0.8 approximately doubled 
the amount of calcium and more than doubled the total solids in the 
leach water. 

Cls did not materially affect the solubility of the grouts made with a 
water-solids ratio of 0.4. Total solubility was increased and calcium 
solubility was decreased for the cls grouts at a water-solids ratio of 
0.8. 

Fly ash did not appear to materially affect solubility of the 0.4 ratio 
grouts. Solubility of the 0.8 ratio grouts seemed to be increased by use 
of fly ash. 

Behavior of cls in the fly ash grouts was the same as in the neat- 
cement grout. 

Slag reduced the solubility of all constituents. The effect on decreasing 
solubility was more marked at the higher water-solids ratio and was 
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more pronounced for the 2 cement: 1 slag than for the 1 cement: 1 slag 
grout. 

h. Pumicite appreciably increased the overall solubility of the 0.4 ratio 
grouts, but slightly decreased the solubility of the 0.8 ratio grouts. 
Calcium solubility was appreciably increased in the 0.4 ratio grouts 
and only slightly increased in the 0.8 ratio grouts. 

i. The opaline shale materially reduced the solubility of all constituents. 
Total solids were reduced by about 50 per cent and calcium by about 

40 per cent. 


SUMMARY OF RESULTS 


The reader should keep in mind in evaluating the data that the investigation 
described comprises only a single set of tests per condition. 


Factors Influencing Penetration of Grouts 


Surface Texture 


It was found that the surface condition of the specimens, that is, whether 
the surfaces to be groutec were smooth or roughened, had a pronounced ef- 
fect on the thickness of the grout that could be forced through the cracks. The 
smoother surfaces permitted the use of considerably lower water-cement 

ratios (thicker grouts) than did the roughened surfaces. 


Water-Cement Ratio 


The water-cement ratio influenced the thickness of crack that could be 
penetrated at a given pumping pressure. Neat grout with a water-cement ratio 
of 0.43 would penetrate the 0.03-in. crack at 25 psi, but the water-cement 
ratio had to be increased to 2.67 before neat grout would penetrate the 0.01- 
in. crack at the same pressure, 


Pumping Pressure 


An approximately straight-line pressure drop occurred along the speci- 
mens while the 0.02- and 0.03-in. cracks were being grouted. The pressure 
was either 25, 50, or 100 psi where the grout entered, and was 0 psi where it 
emerged with intermediate pressures measured along the length of the speci- 
men. The pressure drop was steeper for the first 12 in. of specimen length 
than for the remaining 36 in. when the 0.01-in. crack was being grouted. The 
steepness of the gradient for the first 12 in. increased as the grout thickened 
and the flow decreased. 

Increased pumping pressures caused directly proportionate increases in 
grout flow. At identical water-cement ratios, pressures, and crack thick- 
nesses the neat-cement grout had least flow, cement fly ash grout had more 
flow, and neat grout plus intrusion aid had still more flow. However, the flow 
was not greatest, as might be expected, with cement fly ash grout plus in- 
trusion aid, Data on the cement fly ash intrusion aid grout were anomalous 
and inconclusive. 

When the crack width was 0.02 or 0.03 in. the minimum water-cement ratio 
(thickness of grout) the crack would accept was little affected by the grouting 
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pressure, however, the thickness of the grout that would penetrate the 0.01-in. 
crack depended somewhat on pressure. At 25 psi, the thickest grout that 
would penetrate the 0.01-in. crack had a water-cement ratio of 2.67. When 
the pressure was increased to 50 psi, grout with a water-cement ratio of 1.33 
< would penetrate, but a further increase in pumping pressure to 100 psi did 
not permit use of still thicker grout. 


Grain Size 


The use of a No. 50 sieve cloth for straining the grout seemed beneficial 
and probably the use of a No. 60, or even a No. 100 cloth, if vibrated, would 

q be practical in removing oversized particles which cause trouble in grouting 
the tightest seams. The fly ash contained grains of material larger than 0.01 
in. and grouts containing fly ash would not penetrate the 0.01-in. fissures. 
Based on the maximum grain size of the cement, 0.006-0.01 in., and the thin- 
nest cracks grouted, 0.01 in., the ratio of crack thickness to grain size should 
not be less than 1.7, and probably safer ratios would be 3.0 or more. 

The use of fly ash appears limited to grouting larger cracks unless a sup- 
ply of fine material can be obtained, or unless fly ash can be processed, 
possibly by air separation, to remove the larger granules. No attempt was 
made to grout 0.01-in. cracks with grout containing ground slag, pumicite, 
or opaline shale. Difficulty would doubtless have been encountered with the 
slag and shale whose grain size exceeded 0.01 in. 

The use of more finely ground cements, such as high-early strength 


cement, or ordinary cement processed through an air separator might be 
feasible. 


Chemical Fluidifiers 


The use of intrusion aid appeared to be a help in reducing the water ratio 
of the grout that could be pumped through cracks of 0.03-in. thickness with 
roughened surface, and 0.03- and 0.02-in. thickness with smooth surfaces, It 
is believed that the aid should also have promoted use of lower water-cement 
ratio grout for the roughened 0.02-in. crack but did not in these tests because 
of surface factors in the specimens. Intrusion aid was of no help in promot- 
ing penetration of 0.01-in. cracks at any water ratio tried because of the 
limits of penetration placed upon the grout containing fly ash by the coarse- 
ness of the fly ash. Penetration of 0.01-in. cracks with neat grout containing 
intrusion aid was not tried. The use of intrusion aid appeared to cause a 
small reduction in bleeding and a slight expansion in the grouts. 

Intrusion aid and cls appeared to increase the fluidity of the neat grouts so 

" that a seam of given thickness would pass a greater amount of such grout than 

of plain grout having the same water-cement ratio; conversely, grout with 

fluidifiers, at a lower water-cement ratio and at a given pressure, would 
penetrate a given crack better than grouts without fluidifiers. 

Intrusion aid seemed to be more effective with neat grouts than with grouts 
containing fly ash. The cement fly ash grout with intrusion aid did not pene- 
trate either the 0.02- or the 0,03-in. seams significantly better than the same 
grouts without aid. 


Mineral Fines 


Fly ash appeared to increase the fluidity of the grouts in which it was used. 
The ground slag also appeared to increase fluidity, however this was not 
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conclusive. The pumicite and opaline shale both greatly stiffened the grouts 
in which they were used; however, at comparable consistency with neat grout, 
the pumicite and opaline shale grouts pumped better than the neat grout. At 
equal consistency slag grout pumped about the same as neat grout, but fly ash 
grout did not pump as well. 


Factors Affecting Quality of Grout Films 


Water-Cement Ratio 


The grout films having water-cement ratios of 0.5 or less appeared hard 
and filled the fissures completely with little evident bleeding. However, the 
bond of top to bottom half of the grouted specimen was not considered good 
except when water-cement ratios were lower than 0.4. 

Bleeding occurs in grouts having water ratios in excess of 0.50, and be- 
comes more severe the higher the unit water content. A narrow band of dense 
hard grout was always found next to the shim beneath which water had escaped. 
This band of dense grout was always narrow, never exceeding 3/4 in. in width, 
no matter what the original water content. It has been hypothesized that the 
original water content of a grout is not very important as the excess water 
is squeezed out into the pores of the rock and extremely fine seams (too fine 
to accept grout) leaving a dense hard material filling the cavity. It was not 
found possible to squeeze this water out with the 100-psi grouting pressure 
used, in the 10 to 15 minutes during which pressure was maintained, except 
in a narrow band around the edge of the specimen, even with the outer edges 
of the specimens open to the air at 0 (relative) pressure. However, it should 
be borne in mind that the results under discussion were obtained by grouting 
the space between the two concrete slabs and might not give a true picture of 
what occurs in grouting a foundation hole. 

Setting time increased with water content. Neat grout with a water content 
as high as 2.06 required ten days to set in the setting-time tests. Thin grout 
containing intrusion aid (W/C 4.32) did not set in 14 days. These thin grouts 
can be shown to be incapable theoretically of developing any compressive 
strength when the water-cement ratio of the grout is so high that the space 
between cement grains, which in normal pastes and concretes is of capillary 
proportions and subsequently becomes filled with gel, is so great that it is not 
generally possible for the gel developed by the hydration of the individual 
particles or groups of particles to grow out through the space and make suf- 
ficient contact with that developing from other particles or groups of parti- 
cles.‘3) The higher water-solids ratio grouts were also much less resistant 
to leaching than those with lower ratios. 


Chemical Fluidifiers 


Bonding of the hardened grout films to the top slab was poor (25 per cent 
or less of film adhering to top) for all grouts except those containing fly ash 
with cls where the bond was considered good (40 per cent or more adhering 
to top). Bond strength as measured by shear tests on neat and neat-cement 
plus intrusion aid grouts was relatively low although the use of intrusion aid 
appeared somewhat to increase bonding between top and bottom halves of the 
specimen in neat grout. Cls did not appear to benefit the bond condition of 
specimens intruded with neat grout. 
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Both intrusion aid and cls prolonged setting time and reduced bleeding. 
The solids in the fly ash grouts tended to clump or agglomerate. This tenden- 
cy appeared to be reduced but not eliminated by the fluidifier. Intrusion aid 
was not used in the tests for solubility of grout films, but cls had little effect 
on the leaching characteristics of the grouts tested. 


Mineral Fines 


Hardened grout films containing fly ash showed evidence of channeling 
action and apparent agglomeration. This tendency was not noted when slag, 
pumicite, and opaline shale were used. The use of ground slag, pumicite, and 
opaline shale greatly improved the apparent quality of the grout film, virtually 
eliminating evidence of bleeding, channeling, and agglomeration. Bond be- 
tween top and bottom halves of the specimens was not improved by use of the 
several mineral fines. 

The use of fly ash had little apparent effect on the solubility of the 0.4 
water-ratio grouts, but appeared to increase the solubility of the 0.8 water- 
ratio grout. Slag reduced the leaching in all cases with the most pronounced 
effect in the higher water-ratio grouts. Pumicite increased the overall solu- 
bility of the 0.4 ratio grout but slightly decreased the overall solubility of the 
0.8 ratio grout. The opaline shale greatly reduced the overall solubility of 
both the low- and high- water-solids ratio grouts in which they were used, 


Curing 


Grout films cured for 90 days appeared more resistant to leaching than 
those cured 7 days. 


Measurement of Consistency 


The measurement of consistency was readily accomplished in the labora- 
tory with the torsion viscosimeter. Changes in consistency can also be de- 
tected and controlled readily by means of unit weight measurements. The 
flow cone was not too sensitive to small changes in viscosity, although a cone 
with an orifice restricted enough to require more time to discharge should 
prove suitable for field use. The torsion viscosimeter, because it is delicate 
and requires a quiet, level area for operation, would not be as suitable for 
field as for laboratory use. 


CONCLUSIONS 


1. The surface texture of a fissure has a distinct influence on the thick- 
ness of grout that can be used to fill it. The smoother the surface the lower 
the water-cement ratio of the grout that can penetrate the crack. 

2. The maximum grain size of the solids in the grout determines the mini- 
mum crack width that can be grouted. The ratio of crack width to grain size 
should probably be three or more. 

3. The use of certain chemical admixtures increases the fluidity of grouts, 
to some extent, thereby promoting penetration of a given crack width with 
grouts of a slightly lower water-cement ratio than could be used successfully 
without them. Bleeding can be reduced somewhat and setting time increased 
by use of these materials. 
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4, It was found impracticable to squeeze the excess water from a thin 
grout, thereby leaving a dense hard filler in the cavity, at the pressures and 
with the techniques used in this program. 

5. Bleeding largely prevents bonding of the grouting material to the upper 
surface of the fissure. 

6. The use of finely ground mineral admixtures such as granulated blast- 
furnace slag, pumicite, and opaline shale can reduce the bleeding of a grout 
and greatly improve the continuity and appearance of the hardened grout film. 

7. The solubility of a grout film is influenced by the water-cement ratio 
and composition of the film and length of curing. Certain mineral admixtures 
can be used to reduce the amount of leaching that a grout film undergoes. 

8. A straight-line pressure gradient occurs along a fissure being grouted 
only if the crack is of sufficient width. 
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SYNOPSIS 


This article reports on a study of the properties of glacial lake clays. Data 
collected from a number of sites in the Great Lakes area of the United States 
are presented. Despite local variations in subsoil conditions and the large 
distances between the sites, the clay deposits investigated have much in 
common and it is possible to study them as a related family. An analysis was 
made of the principal factors that control the strength characteristics of the 
clays. The shearing resistance and sensitivity were found to be intimately 
related to the structure of the clay as well as the stress history. 


INTRODUCTION 


In North America the activities of continental glaciers led to the formation 
of many glacial lakes during Pleistocene times. Many of them have since 
disappeared or now exist in diminished sizes. A number of soft clay deposits 
extending to great depths are known to exist in areas that had been formerly 
occupied by glacial lakes. Among the most notable are the clays in the Great 
Lakes region such as the deposits ai Chicago and Detroit. These deposits 
consist primarily of a heterogeneous and unstratified silty clay with some 
sand and gravel. The texture of these soils resembles closely that of glacial 
tills. Also present in these areas is the laminated clay and silt known as 
varved clay and generally considered as representative of lacutrine sediments. 

As the majority of the glacial lake clays exhibit low strength accompanied 
by high compressibility, their properties have long been subjects of interest. 
A number of significant publications have described the engineering proper - 
ties of lake clays as well as regional subsoil conditions. This article contains 
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information on the physical properties and strength characteristics of several 
clay deposits in the Great Lakes region. Because of their uniformity and ex- 
tensive depth, these deposits were selected for detailed study. The funda- 
mental factors that govern the important characteristics as sensitivity and 
strength-consolidation relationship of these deposits were analyzed. Also in- 
cluded in this article are some supplementary data from a number of sites. 
These, although lacking in completeness, provide considerable information 
that support the main thesis. 


Description of Clay Deposits 


During the last great ice invasion, the Wisconsin, the continental glacier 
moved into the area from the north and east. At its climax, the glacier 
covered all of Michigan and extended deep into Indiana and Ohio. Subsequent 
recession of the glacier followed the same general directions. The numerous 
retreats and readvances of the glacier resulted in the formation, at various 
times, of large ice-front lakes. The shoreline of these lakes are traced on 
the map in Fig. 1.(1,2) 

The four soil deposits investigated are all in the Great Lakes area. They 
are situated respectively, near the cities of Freemont, Detroit, Saginaw and 
Sault Ste. Marie as indicated in Fig. 1. Fig. 2 is a profile showing the subsoil 
conditions north of Freemont, Ohio. The major units in the subsoil consist of 
a layer of stratified clay and silt overlying a deposit of unstratified silty clay 
with sand and gravel. Both are soft to medium in consistency. Beneath these 
are found one or more old till sheets of very stiff clay with gravel and some 
boulders. Bedrock is encountered at elevations of about 540. Representative 
boring logs describing the physical properties of the soft clays at two sites 
are shown in Figs. 3 and 4. 

Very similar subsoil conditions are found in the Detroit area. Fig. 5 con- 
tains a detailed log and summary of the test results obtained from a site near 
the Industrial Expressway and River Rouge. The data are representative of 
the subsoil conditions in this area. The uppermost unit is a layer of stratified 
clay and silt and it lies over a fairly thick deposit of unstratified clay with 
sand and gravel. In the transition a thin stratum of clay containing many ir- 
regular layers and lenses of sand is revealed in some borings. The thickness 
of this stratum where encountered, does not exceed two feet. At a depth of 
75 ft., the unstratified clay with sand and gravel changes gradually into a 
layer of soft clay with occasional silt laminations. Remnants of old till sheets 
and heavily precompressed lacustrine clay are found between the soft clay 
and bedrock, which is located at about elev. 480. 

At the site near Saginaw, the unstratified silty clay extends down to depth 
of about 51 feet. A layer of stratified clay and silt is found between depths of 
51 and 65 feet and beneath this are layers of sand and still clay till. Data ob- 
tained from a representative boring are shown in Fig. 6. 

Between Mackinac Straits and Sault Ste. Marie, in Northern Michigan is an 
extensive deposit of soft stratified clay and silt up to 60 ft. thick at some lo- 
calities. The presence of hematite gives the clay a distinctive red color. 
Borings were made at two sites near Sault Ste. Marie and typical test results 
from each site are presented in Figs. 7 and 8. 

Detailed study was made of the clay fraction of all the soils encountered 
by means of x-ray diffraction and differential thermal analyses. The results 
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indicated that the principal constituent minerals are quartz, kaolinite, and 
illite, with small percentages of chlorite and vermiculite. The general uni- 
formity in mineralogical contents is well expressed in Fig. 9, in which the 
plasticity indices of the clays are plotted against the clay fractions. All the 
clays belong to the class of inactive clays (3) The index properties of the 
clays are summarized in Table 1. 


Strength Characteristics 


A comprehensive series of laboratory tests was performed to evaluate the 
strength and stress history of the clays. The tests included the consolidation, 
unconfined compression and triaxial tests. The subsequent paragraphs 
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describe the details of the laboratory investigation. A summary of the 
strength properties is given in Table 2. 


Specimens 


All specimens used in the laboratory investigation were obtained by seam- 
less steel tube samplers with a wall thickness of 1/16 inch. From the 
Fremont site only two-inch diameter specimens were obtained. Both two-inch 
and three-inch diameter specimens were used in the investigation of the 
Detroit site. In this case, consolidation tests were performed only on the 
large diameter specimens, whereas, compression tests were performed on 
specimens of both sizes. The specimens from Saginaw were taken with three- 
inch tube samplers and the specimens from Sault Ste. Marie were taken with 
three-inch piston samplers. 


Consolidation Tests 


Some typical results of standard consolidation tests are plotted in Figs. 10 
and 11. From these pressure vs void-ratio curves the preconsolidation 
pressures were determined by Casagrande’s graphical method. Figs. 3 
through 8 show the measured preconsolidation pressures as well as the exist- 
ing overburden pressures at the various sités. The good agreement between 
the measured preconsolidation pressures and the existing overburden indicate 
the most of the soft to medium clay deposits are normally loaded. Only the 
top part of the deposits exhibit substantial precompression and these can be 
attributed to desiccation due to fluctuations in the local ground water table. 


Shear Tests 


The shear strength of the clays was measured by unconfined compression 
tests and undrained triaxial tests. The undrained triaxial tests utilized an 
all-around pressure equal to the value of the existing overburden pressure. 
The deviator stress was applied at a constant rate of strain of 2% strain per 
minute. The same rate of strain was adopted for the unconfined compression 
tests. When three-inch specimens were available, they were trimmed down 

to a diameter of 2.8 inches for compression tests. The 2 inch specimens were 

not trimmed. All specimens have a length-diameter ratio of 2.0. Since speci- 
mens of both sizes were tested for the Detroit site an attempt was made to 
evaluate the effect of the specimen size on strength. Comparison of results 

from five borings revealed no noticeable difference in strength between the 

two types of specimens. 

In Figs. 3 through 8 are plotted the values of compressive strength as de- 
termined by the laboratory tests. The compressive strength was taken as : 
equal to the maximum deviator stress. At all sites there is a general tenden- 
cy for the strength to increase with depth within each individual stratum. 

Near the ground surface, the effect of preconsolidation is reflected in the 

relatively high strengths. However, some exceptions are noteworthy. At 

Detroit, a thin layer of soft clay is encountered at a depth of 75 ft. This is 
designated as clay no. 3 in Fig. 5. A similar clay exists at Saginaw between 

depths of 51 and 65 ft. Both units possess strengths far below those of the 
overlying soils. The sensitivity of the clays varied from 2 to almost 10. The 
sensitive clays include the clays at Sault Ste. Marie, clay no. 3 at Detroit and 

clay no. 2 at Saginaw. All the other clays are classified as having low 

sensitivity. 
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The Ratio 


Because of the high degree of uniformity of the soil deposits investigated, 


it is possible to analyze the relationship between strength and consolidation 
pressure in further detail. 


From the consolidation and shear test data the ratio . was computed. In 


this ratio c denotes the shear strength or one-half the compressive strength 
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Void Ratio 


Pressure - kg /cm* 


Figuee |/ CONSOLIDATION CURVES - SAULT STE. MARIE CLAY 


and p the preconsolidation pressure. For normally loaded clays, the existing 
overburden pressures were used in the computations. The computed values 
are listed in Table 2. 

Besides the clay deposits reported in the preceding paragraphs considera- 
ble data are also available from other sites in this area, especially Chicago. 
The Chicagoclay is an unstratified silty clay with variable percentages of 
sand and gravel and extends to a depth of about 60 ft.(4) G. HL Otto, after a 
thorough analysis of numerous boring data, concluded that the Chicago clay 
consists of a succession of till sheets deposited under water.‘5) Consolidation 
test results presented by J. D. Parsons and R. B. Peck(6) indicate that the 
range of measured preconsolidation pressure agrees well with the existing 
overburden between elevations of -5 and -22, Chicago Datum. The most 
probable preconsolidation pressure as determined by the authors, however, 


/8 
Spec./-7-5 (48) 
an 


ASCE GLACIAL LAKE CLAYS 1732-15 


is slightly greater than the existing overburden. These results are in good 
agreement with those from another site in East Chicago reproduced in Fig. 12. 


When the most probable consolidation pressure is used, the computed £ ratio 
for the Chicago clay is about 0.20. P 

In the Cleveland area the repeated advances and retreats of the glacier ice 
alternated with the formation of ice front lakes. This has resulted in a com- 
plex series of lacustrine clays with or without silt laminations and unstrati- 
fied clays with sand and gravel classified as tills. The Cleveland clay is 
heavily preloaded and no particularly soft materials are found in the area, 
From the data obtained by Peck(7) it is possible to estimate the value of the 


- ratio as somewhere between 0.17 and 0.18 for units Band D. Fig. 13 con- 
tains laboratory test results on clay unit A at a site in the same area. The : 


ratio is about 0.19. 

At Green Bay, Wisconsin, a thick deposit of reddish clay exists to a depth 
of over 100 ft. The deposit consists of a succession of reddish silty clay till 
and some stratified material of lacustrine origin. Data from one boring are 
presented in Fig. 14. Between elevations 530 and 550, the clay is uniform. 
Preconsolidation pressures were determined only at elevations 539 and 540. 


From this, the computed © ratio is 0.19. 
In Fig. 15 the : ratios are plotted against the plasticity indices. The 


straight line in the figure is the relationship obtained by A. W. Skempton(8) 
and L. Bjerrum(9) for estuarine and marine clays. The overall trend is one 


of increasing : ratio with increasing plasticity index. However, the data ob- 


tained in this investigation appear to fall into two distant alignments. The 
first one, consisting of clays with low sensitivity, is slightly above but very 
close to Skempton’s line. The second group falls far below the first one, and 
is composed entirely of sensitive clays. 

The liquidity indices of the clays are plotted against their preconsolidation 
pressures in Fig. 16. Also shown in the figure are curves for the Horten 
clay(10) with a sensitivity of 17 and the Gosport clay(11) with a sensitivity of 
2.4. For a given clay, consolidation under increasing pressures always pro- 
gressively reduces the liquidity index. However, as is obvious in Fig. 16, 
different clays may possess very different liquidity indices at same pre- 
consolidation pressures. As explained by A. W. Skempton(12) this relation- 
ship depends primarily upon the microstructure of the clay particles. Since 


‘ the data in Fig. 13 show that low ~ ratio and high sensitivity are limited to 


clays with a high liquidity index, it follows that these properties must be de- 
termined by the micro structure of the clay. 


Stress-Strain Relationship 


The larger number of unconfined compression tests also yielded consider- 
able information on the stress-strain characteristics of the clays. This re- 
lationship may be classified into two distinctive types as illustrated in Fig. 17. 
The plot for specimen 5-13-4 Detroit exhibits a brittle material with a rather 
high value of initial modulus. The stress rapidly decreases after attaining 
a peak value at low strain. In contrast to this, the plot for specimen 5-5-4 
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Fig. 16 RELATIONSHIP BETW. PRECONS. PrEssVeE &L/ 


Detroit has the characteristic plastic nature. Its maximum stress is attained 
at a large strain of 35% and there is no significant reduction in stress with 
further strain, Photographs of the failed specimens are shown in Fig. 18. 
The strain at which maximum strength is attained with the various clays are 
listed in Table 2 as Ep: The brittle soils are all sensitive clays and all 
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(a) Brittle Failure 
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evidence indicates that the stress-strain characteristics are governed by the 


same factors that govern the © ratio and sensitivity. 


Fabric of Clays 


The investigation on strength has revealed that there are two groups of 
clays with very different strength characteristics in spite of their close simi- 
larity with respect to activity and mineralogical contents. It is further indi- 
cated that the major differences may be attributed to the micro-structure of 
the clays. The attention was, therefore, turned to the fabric of the clays. The 
fabric of a clay is usually defined as the arrangement of the clay flakes inside 
a given specimen or deposit. This arrangement may be either parallel or 
random depending upon the depositional environments. 


Optical Studies 


The optical properties of the clay minerals provide some means to evaluate 
the orientation of the clay flakes. Although most clay minerals are biaxial 
crystals, the difference in the indices of refraction are known to be small; and 
for the present study they are treated as uniaxial crystals with optical axes 
represented by aj, ag, and c, with the c axis perpendicular to the flat face of 
the particle as shown in Fig. 19a. 

When flake-shaped clay-mineral particles settle from a deflocculated sus- 
pension, they come to rest with flakes one on top of another so that their basal 
plane surfaces are more or less parallel. Mitchel has found this to be true 
also for fresh water sediments.(13) In such a clay, the c axes of all the flakes 
are oriented in the vertical direction although the aj and ag axes of the flakes 
may not be parallel (Fig. 19b). The orientation can, therefore, be studied by 
examination of thin sections under a petrographic microscope. A vertical 
section of an oriented clay will exhibit four alternate stages of illumination 
and extinction when placed between crossed nicols of the microscope since 
light passing through such a section is polarized to vibrate in two perpendicu- 
lar planes. A horizontal section will show no extinction or illumination as the 
ay and ag axes of the individual flakes lie in all directions (Fig. 19c). On the 
other hand, in a group of randomly oriented particles, none of the optical axes 
are alined in any given direction and illumination and extinction are non- 
existent in both vertical and horizontal sections. The predominance of the 
stages of illumination and extinction revealed in a vertical section, therefore, 
can be used as an index of the degree of orientation of the clay particles. 

In this study the intensities of light transmitted through the microscope at 
the stages of illumination and extinction were measured by a phototube. The 
ratio of the light intensity at illumination to that at extinction was determined 
readily. This ratio is designated the orientation factor and reflects the ar- 
rangement of clay particles of the specimen. This ratio should be equal to 
infinity when a specimen exhibits complete extinction and is equal to one if 
the stages of extinction and illumination are non-existent. This sections 30 
microns thick were used. The thin sections were prepared from clay speci- 
mens impregnated with a polyethylene glycol (Carbowax 6000) following the 
procedure described by J. K. Mitchell. 13) At least two thin sections were 
made from each clay and light readings were taken through the microscope 
with a 10x eyepiece and objectives of 1.0x, 10x and 45x, The values of the 
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readings are listed in Table 3. However, when the oriented areas were very 
small and when the directions of orientation were not the same in all areas, 
it was possible only to estimate the degree of orientation by visual exami- 

nation. 


BS 


Significance of Fabric 


On the basis of the optical studies it is possible to further classify the 
clays according to their fabric. The first conclusion that may be drawn upon 
an examination of the data in Table 3 is that the till-like, unstratified silty 
clays all show random orientation over a large area. (Orientation factor = 1). 
This, of course, is understandable considering the large percentages of sand 
and silt in the soil. When magnified 450 times, it was found that the clay is 
actually well oriented around the silt and clay particles as shown in Fig. 20. 

The fabric of the stratified lacustrine clays varied from well oriented (ori- 
entation factor from 2 to 3) to almost random. Fig. 21 shows micrographs of 
the illumination and extinction stages in a well oriented specimen and Fig. 22 
micrograph of a specimen with little orientation. Both Detroit clay no. 1 and 
Fremont clay no. 1 exhibit a high degree of parallel arrangement whereas the 
Sault Ste. Marie clays 1 and 2 and Saginaw clay no, 2 show very little orien- 
tation. Since the percentages of silt and sand is very low in these clays, this 
lack of orientation can only result from a random structure. The findings re- 
garding the particle orientation further substantiate the concept of micro- 
structure discussed previously. A flocculent or honeycomb structure is 
necessarily accompanied by random particle orientation, It is, therefore, this 
structure that gives the clays their high liquidity index, high sensitivity, and 


low : ratio. The clays with pronounced parallel particle arrangement and the 
tills have low sensitivity and their £ ratios are close to values reported in the 
works of A. W. Skempton(8) and L. Bjerrum. (9) 


Conditions of Deposition 


To further check the micro-structure of the clays, the concentration of 
soluble salts in the porewater and the exchangeable cations were determined 
for the clays. It was reasoned that the flocculent structure of some clays may 
be the result of flocculation during deposition. This conceivably could have 
been brought about by deposition in waters of high salt concentration or by 
the very nature of the exchangeable cations. (14) 

The results of soluble salt measurements and chemical analyses are given 
in Table 4. The salt concentration of lacustrine clays are very low and do 

€ not differ materially between clays with widely different structures. Further- 
more, there is no significant difference in the exchangeable cations. The 
difference in micro-structure exhibited by the clays, therefore, cannot be 

explained by the results of chemical analyses. 


CONCLUSIONS 
In summary it may be said that the soft clays investigated are either la- 
custrine clays or tills deposited by glaciers under water. Consolidation of 
these soils is accomplished largely by the weight of the overburden and most 
of the clays can be considered as normally loaded. Despite the overall 
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(b) Extinction 


Micrographs of a specimen with poorly oriented clay particles 


Sault Ste. Marie Clay No. 1 (Mag. 100 X) 


(a) Illumination 
Figure 22. 
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similarity with respect to mode of deposition and mineralogical content and 
micro-structure of the clays was found to vary widely. In the case of la- 
custrine clays, the structure may be honeycombed and flocculent or the ar- 
rangement of the flake-shaped clay particles may be more or less parallel. 
The structure is accurately reflected by the relationship between the liquidity 
index and preconsolidation pressure. Important strength characteristics as 


sensitivity and the - ratio were found to be dependent upon the micro-structure 


of the clays. 
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Table 2. Strength Characteristics of Clay Deposits 


Location Unit P(Kg/cm*) c/P 


August, 1958 


F remont 1.00 
1.60 


Detroit 9.00 
1.70 
2.65 


Saginaw 


Sault Ste. Marie 
Chicago!) 


Cleveland’ 


Green Bay 


p - preconsolidation pressure 
c - 1/2 compressive strength 
sensitivity 

€ - strain at max. stress 
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